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19  ABSTRACT  (Continued) 


the  field  inside  a  continuous  enclosure  is  a  monotonically 
decreasing  function  of  frequency.  At  frequencies  above  cutoff, 
the  fields  in  both  slotted  and  continuous  enclosures  show  very 
complex  behavior  associated  with  cavity  resonances. 

The  frequency  dcmain  series  for  the  internal  magnetic  fields 
were  transformed  into  time  danain  expressions  which  were  used  to 
plot  the  transient  response  of  the  internal  magnetic  fields  at 
various  locations  when  one  wall  or  one  slot  is  exposed  to  time- 
varying  external  fields.  It  was  found  that  the  properties  of  the 
field  inside  a  continuous  enclosure  are  determined  to  a  great 
extent  by  the  characteristic  diffusion  time  of  the  enclosure  wall 
and  the  duration  of  the  incident  field.  The  former  determines  the 
rastcst  rise  time  of  the  internal  field  and  the  latter  determines 
its  largest  value.  The  transient  field  inside  an  enclosure  with  an 
air-filled  slot  consists  of  the  sum  of  two  parts:  the  stationary 
field  and  the  propagating  field.  The  stationary  field  dominates 
close  to  the  slot  where  it  forms  a  reduced  replica  of  the  external 
field.  Its  peak  field  decreases  rapidly  with  distance  from  the 
slot  so  that  at  most  locations  it  is  negligible  compared  to  the 
propagating  field.  The  latter  travels  at  approximately  the  speed 
of  light  and  undergoes  repeated  reflections  from  the  walls  of  the 
enclosure.  The  peak  value  of  the  stationary  field  depends 
principally  on  the  length  of  the  slot,  while  the  peak  of  the 
propagating  field  depends  on  the  rise  time  of  the  external  field. 
That  is,  longer  slots  produce  larger  stationary  fields  and  faster 
rise  times  produce  larger  propagating  fields. 

Calculations  based  on  the  frequency  domain  expressions  were 
used  to  compare  the  theoretical  results  with  two  sets  of  experi¬ 
mental  data,  and  satisfactory  agreement  was  found. 
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INTRODUCTION 


1  . 


An  electromagnetic  shield  is  a  structure  fabricated  from 
one  or  more  sheets  of  electrically  conducting  material,  usually 
copper,  steel,  or  aluminum,  that  is  used  to  reduce  the  intensity 
of  electric  and  magnetic  fields  entering  a  designated  region  (the 
shielded  volume)  from  external  sources  such  as  antennas,  trans¬ 
mission  lines,  lightning,  and  electromagnetic  pulse  ( EMP ) .  It 
does  this  by  forming  a  barrier  between  the  shielded  volume  and 
the  source  that  reflects  and  absorbs  most  of  the  electromagnetic 
field  before  it  can  enter.  Shields  can  be  simple  open  structures 
such  as  a  single  flat  sheet  of  steel  placed  between  a  source  and 
the  shielded  volume.  Or,  they  can  be  much  more  complicated  closed 
form  structures  that  approximate  continuous  metallic  shells.  The 
latter  are  usually  circular  cylinders  or  rectangular  parallelepi- 
pipeds  (boxes).  The  volume  enclosed  by  such  a  shield  may,  or  may 
not,  coincide  with  the  shielded  volume.  It  does  if  the  source  of 
interest  is  located  outside  of  the  enclosure.  However,  if  the 
source  is  inside,  then  the  shielded  volume  is  the  (unlimited) 
region  outside  the  enclosure.  On  the  other  hand,  if  there  are 
sources  of  interest  both  inside  and  outside,  then  both  regions 
are  shielded  volumes  with  respect  to  one  source  or  the  other. 
That  is,  a  closed  shield  may  be  used  to  exclude  designated  fields 
from  a  limited  region,  it  may  confine  fields  to  a  limited  region, 
or  it  may  separate  the  effects  of  two  sources  by  excluding  the 


fields  of  one  source  from  the  region  while  confining  the  fields 
of  the  other  source. 

Closed-form  shields  are  prefered  to  open  shields  for  most 
applications  because  of  their  generally  greater  efficiency  as 
measured  by  the  reduction  in  field  intensity  at  a  given  location 
when  that  location  is  included  in  a  shielded  volume.  The  reason 
for  this  is  simply  that  they  usually  have  smaller  open  (non- 
metallic)  areas  through  which  fields  can  reach  the  shielded  vol¬ 
ume  without  incurring  the  large  reflection  and  absorption  losses 
associated  with  direct  penetration  of  the  metal  barrier.  Open 
shields,  by  their  very  openness,  usually  provide  relatively  large 
areas  where  incoming  fields  can  bypass  the  metal  barrier.  For  ex¬ 
ample,  a  shield  consisting  of  a  flat,  rectangular,  steel  sheet  of 
modest  thickness  exacts  huge  losses  on  the  part  of  the  incoming 
field  that  takes  a  straight  line  all-metal  path  to  the  shielded 
volume.  If  this  were  the  only  path  to  the  shielded  volume,  its 
efficiency  would  be  very  high.  However,  this  structure  also  pro¬ 
vides  longer  but  less  lossy  paths  through  the  (infinite)  area  at 
and  beyond  the  perimeter  of  the  sheet  by  which  part  of  the  field 
can  reach  the  shielded  volume  through  the  process  of  diffraction. 
This  diffracted  field  usually  dominates  the  directly  transmitted 
field  and,  thus,  severely  limits  the  efficiency  of  the  shield. 
Without  changing  the  frequency  of  the  source  field,  the  only  way 
the  efficiency  can  be  improved  is  by  increasing  the  size  (area) 
of  the  sheet-  an  approach  that  is  often  impractical. 


A  simple  closed-form  shield  can  be  constructed  from  the 
rectangular  steel  sheet  by  bending  it  into  a  circular  cylinder 
and  joining  the  adjacent  edges  with  a  continuous  weld.  Since  the 
diffraction  field  can  now  avoid  all-metal  paths  only  by  entering 
or  leaving  the  enclosed  volume  through  the  finite  (circular) 
areas  at  the  ends  of  the  cylinder,  we  can  expect  the  efficiency 
of  this  shield  to  be  significantly  greater  than  that  of  a  flat 
sheet  with  the  same  composition,  thickness,  and  surface  area. 
Much  greater  efficiency  can  be  achieved  by  welding  steel  sheets 
to  both  ends.  These  minimize  diffraction  effects  and  insure  that 
fields  enter  or  leave  the  enclosure  by  all-metal  paths.  In  this 
way,  a  closed-form  shield  can  be  transformed  into  a  continuous 
metal  shell. 

By  a  similar  process  in  which  six  rectangular  steel  sheets 
are  welded  along  their  edges,  a  continuous  metal  shell  can  be 
constructed  in  the  form  of  a  rectangular  parallelepiped.  This 
enclosure,  like  the  cylindrical  enclosure,  constitutes  an  ideal 
electromagnetic  shield  in  the  sense  that  its  efficiency  cannot 
be  significantly  improved  without  increasing  the  thickness  of  the 
shell  or  changing  its  composition.  It  does,  in  fact,  realize  the 
full  shielding  potential  of  a  given  thickness  of  steel  in  this 
form.  This  structure  is  also  ideal  in  the  sense  that  it  repre¬ 
sents  a  limit  that  can  only  be  approached  by  practical  shields. 


In  practice,  perfect  metallic  continuity  must  always  be 
sacrificed  for  access  areas  in  the  form  of  doorways,  hatches,  and 
other  entry  points  for  people,  equipment,  or  supplies.  These 
openings  or  discontinuities  inevitably  reduce  the  efficiency  of  a 
practical  shield  below  that  of  the  ideal  by  allowing  fields  to 
reach  the  interior  by  nonmetallic  paths  just  as  in  the  case  of  an 
open  shield.  Thus,  practical  closed-form  shields  must  be  open  to 
some  extent.  To  minimize  their  effect,  large  openings  are  usually 
provided  with  metal  closures  (doors  or  panels)  that  are  intended 
(when  in  place)  to  establish  continuity  by  maintaining  metal  to 
metal  contact  at  their  boundaries  with  the  rest  of  the  shield. 
Unfortunately,  perfect  contact  between  two  pieces  of  metal  cannot 
usually  be  maintained  except  by  welding  or  some  other  type  of 
permanent  connection  that  would  defeat  the  purpose  of  the  open¬ 
ing.  Consequently,  closures  are  almost  never  completely  success¬ 
ful  in  achieving  continuity,  although  some  may  approach  it.  In 
general,  the  effect  of  the  closure  is  to  replace  an  original 
large  opening  by  one  or  more  smaller  ones  in  the  form  of  gaps 
arranged  along  the  seams  where  it  meets  the  body  of  the  shield. 
These  openings  are  frequently  sufficient  to  reduce  the  efficiency 
of  the  structure  well  below  that  of  a  comparable  ideal  shield, 
and,  if  numerous  enough  and  large  enough,  are  capable  of  compro¬ 
mising  its  performance.  Consequently,  the  treatment  of  seam  dis¬ 
continuities  is  crucial  to  both  the  theory  and  practice  of  elec¬ 
tromagnetic  shielding 


This  report  presents  a  theoretical  investigation  of  the 
rectangular  enclosure  as  an  ideal  shield  and  as  a  shield  with  one 
or  more  seam  discontinuities  in  the  form  of  narrow  rectangular 
slots.  The  investigation  uses  a  recently  developed  shielding 
theory  for  generalized  structures  composed  of  flat  metal  sheets 
to  obtain  general  expressions  for  the  electric  and  magnetic 
fields  at  points  inside  continuous  and  discontinuous  rectangular 
enclosures  when  an  external  surface  is  exposed  to  fields  from  an 
arbitrary,  time  harmonic,  electromagnetic  source.  These  general 
expressions  are  reduced  to  explicit  expressions  when  the  source 
field  is  spatially  uniform  over  the  outside  surface  of  the  enclo¬ 
sure  and  for  the  less  restrictive  case  where  the  field  is  uniform 
over  the  discontinuity.  The  latter  are  used  to  calculate  and  plot 
frequency  domain  fields  at  selected  points  inside  both  continuous 
and  discontinuous  enclosures.  Comparison  of  the  calculations 
shows  the  large  effect  of  even  small  discontinuities  on  the 
shielding  effectiveness  (efficiency)  of  the  enclosure.  Analogous 
results  are  obtained  in  the  time  domain  by  computing  and  plotting 
inverse  LaPlace  transforms  of  the  frequency  domain  expressions 
for  several  types  of  transient  source  fields.  In  the  case  of  the 
discontinuous  enclosure  these  calculations  reveal  a  very  compli¬ 
cated  time  history  for  the  internal  fields  due  to  multiple  inter¬ 
nal  reflections. 

The  following  section  reviews  the  general  shielding  theory 
based  on  impedance  boundary  conditions  and  summarizes  its  princi¬ 
pal  results.  Section  3  uses  these  results  to  construct  the  gen- 


lJVmV-V > -VVVVHVV..N  -V.  •  .V.V-  A  A.  -  A  A  A  A  A  A  A  A  A  A  A  .  •  A  A  .  •  A  A  A  A  A  A  -A  A  c.  A  A  A 


'  «  %  "  -  *  .  •  t*  ..■  ■/  -/  V*  •/  •. 


eral  expressions  for  internal  electric  and  magnetic  fields  as 

doubly  infinite  series  summed  over  the  TE  modes  for  a  rectangu- 

nm 

lar  waveguide.  Fourier  coefficients  appearing  in  the  series  are 
evaluated  in  section  4  for  the  cases  of  interest,  and,  in  section 
5,  the  frequency  domain  series  are  transformed  to  the  time  domain 
using  term-by-term  inversion  .  In  section  6,  calculations  are 
carried  out  using  some  of  the  expressions  obtained  in  the  preced¬ 
ing  sections,  and  these  are  compared  to  measurements  performed  by 
two  groups  of  experimental  investigators. 


2. 


A  SHIELDING  THEORY  BASED  ON  IMPEDANCE  BOUNDARY  CONDITIONS 


In  1985,  Monroe*  described  a  general  theory  of  electromag¬ 
netic  shielding  that  improved  on  and  extended  the  classic  trans- 

2 

mission  theory  of  Schelkunoff.  The  new  theory  improved  on  the 

old  one  by  replacing  its  transmission  line  model  and  eliminating 

certain  unnecessary  assumptions.  The  transmission  line  model  was 

3  4 

replaced  by  an  Impedance  Boundary  Condition  (IBC)  ’  which  had 
not  previously  been  used  explicitly  to  solve  shielding  problems. 
This  approximate  boundary  condition  performs  the  same  function  as 
the  transmission  line  model  in  that  it  allows  one  to  represent 
the  field  passing  through  a  wall  of  the  shield  as  a  plane  wave 
traveling  perpendicular  to  that  wall.  It  also  provides  a  way  to 
calculate  this  field  in  terms  of  the  source  field  incident  on  the 
outside  of  the  shield.  Its  advantage  over  the  transmission  line 
model  is  that  it  can  be  applied  to  a  very  much  larger  class  of 

*R.  L.  Monroe,  A  Theory  of  Electromagnetic  Shielding  with 
Applications  to  MIL-STD  285,  IEEE-299,  and  BMP  Simulation, 

Harry  Diamond  Laboratories,  HDL-CR-85-052- 1 ,  Adelphi,  MD 
(February  1985). 

2 

S.  A.  Schelkunoff,  Electromagnetic  Waves,  Van  Nostrand,  Prince¬ 
ton,  NJ  (1943). 

^T .  B.  A.  Senior,  Appl .  Sci.  Res.,  8  (B)  (1960),  418. 

4 

T.  B.  A.  Senior,  IEEE  Trans,  on  Antennas  Propagt.,  AP-29,  No.  2 


sources  and  shields,  including  shields  that  have  been  rendered 
discontinuous  by  one  or  more  slot-like  apertures.  In  fact,  the 
restrictions  on  applying  IBC’s  are  so  weak  that  there  is  virtual¬ 
ly  no  electromagnetic  shield  to  which  an  IBC  could  not  be  appli¬ 
ed.  Thus,  with  the  aid  of  IBC’s,  one  can  develop  a  theory  of 
shielding  for  a  wide  range  of  sources  and  shields  that  fall  out¬ 
side  the  scope  of  Schelkunof f ’ s  theory. 


2.1  Impedance  Boundary  Conditions 


An  IBC  is  a  relationship  between  an  impedance  func¬ 
tion  and  an  electromagnetic  field  at  the  interface  between  two 
electrically  distinct  media  where  the  impedance  function  charac¬ 
terizes  one  medium  and  the  electromagnetic  field  is  defined  in 
the  other  medium.  In  its  most  frequently  applied  form,  the  IBC 
relates  tangential  field  components  at  the  interface  to  the  im¬ 
pedance  looking  into  one  of  the  media.  With  the  two  media  labeled 
Ml  and  M2  as  shown  in  figure  1,  this  condition  can  be  written  in 
vector  form  as  follows: 


n  x  (n  x  El)  =  -  rjZ  (n  x  HI) 


(2.1) 


where  El  and  HI  are  t’i®  electric  and  magnetic  fields  in  Ml ,  n  is 
the  unit  vector  normal  to  the  interface  pointing  outward  from  M2, 
q 2  is  the  impedance  looking  into  M2,  and  it  is  understood  that 
(2.1)  applies  only  at  the  interface.  Equation  (2.1)  is  an  approx¬ 
imation  that  cannot  be  used  to  replace  standard  boundary  condi¬ 
tions  in  the  general  case.  However,  in  many  cases  of  interest,  it 
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x  El  and  n  x  HI  at  the  boundary  between  Ml  and  a  second 


medium  M2  act  as  a  primary  source  for  the  field 


has  been  shown  that  (2.1)  is  a  valid  approximation,  and,  in  these 
cases,  the  IBC  can  be  used  to  simplify  the  problem  of  determining 
the  fields  in  Ml  and  M2. 

This  simplification  derives  from  the  fact  that  (2.1)  de¬ 
couples  the  fields  in  Ml  from  those  in  M2  in  a  way  that  does  not 
introduce  spatial  derivatives  of  the  fields  at  the  interface. 

This  means  that  El  and  HI  can  be  computed  independently  of  the 

fields  (E2,H2)  in  M2  and  that  both  sets  of  fields  can  be  obtained 
by  applying  standard  techniques  to  Maxwell’s  equations.  One  first 

solves  Maxwell’s  Equations  for  El  and  HI  using  (2.1)  to  replace 

M2  and  then  solves  for  E2  and  H2  in  M2  using  El  and  HI  at  the 
interface  to  replace  Ml.  Since  this  two-step  process  will  usually 

be  much  easier  than  solving  Maxwell’s  equations  directly  for  El, 

HI,  E2,  and  H2 ,  the  utility  of  (2.1)  is  obvious.  Moreover,  if  one 

is  interested  only  in  the  fields  in  Ml,  then  E2  and  H2  need  not 
be  computed  at  all.  Conversely,  if  one  is  primarily  interested  in 

E2  and  H2 ,  then  it  is  only  necessary  to  solve  for  El  and  HI  at 
the  interface  in  order  to  determine  the  fields  throughout  M2.  The 
latter  describes  the  usual  situation  in  shielding  problems  where 
the  interior  of  M2  can  be  identified  with  the  shield,  Ml  is  the 
region  containing  one,  or  more,  electromagnetic  sources,  and  only 
the  fields  transmitted  into  the  shielded  volume  are  of  interest. 


To  take  advantage  of  the  IBC,  it  is  necessary  to  establish 
the  validity  of  (2.1)  at  the  interface  that  defines  the  problem 

of  interest.  In  general,  this  requires  one  to  show  that  E2  and  H2 
propagate  into  M2  along  n  as  approximate  plane  waves.  One  way  to 

do  this  is  to  show  that  the  variation  of  E2  and  H2  along  n  is 

much  larger  than  the  variation  of  El  and  HI  at  the  interface  in 
directions  transverse  to  n.  Specifically,  one  can  show  that  the 

normal  derivatives  of  E2  anf  H2  are  much  larger  in  magnitude  than 

the  transverse  derivatives  of  El  and  HI  at  the  interface.  Since 
electromagnetic  shields  are  intended  to  produce  just  such  varia¬ 
tions,  it  is  not  suprising  that  most  of  these  structures  will  be 
found  to  satisfy  IBC’s.  Indeed,  it  can  be  argued  that  no  struc¬ 
ture  can  be  an  effective  shield  unless  it  does  satify  an  IBC  at 
all  points  on  its  surface. 

The  validity  of  (2.1)  has  been  demonstrated  under  relative¬ 
ly  weak  restrictions  in  the  case  of  a  planar  interface  separating 
two  homogeneous  half-spaces  as  shown  in  figure  2,  where  Ml  is 
free  space  and  M2  consists  of  a  material  with  complex  permittiviy 
t2  and  permeability  fj2 .  With  this  arrangement,  equation  (2.1)  can 
be  written  in  scalar  form  as  follows: 

El  =  -  172  HI  ,  El  =  rj2  HI  at  z  =  0,  (2.2) 

x  y  y  * 

where 

1/2 


qZ  =  (fj  2/t2) 


(2.3) 


and  the  fields  are  referred  to  a  rectangular  coordinate  system 
with  its  origin  at  the  interface  and  its  z  axis  directed  out  of 

A 

M2  (n  =  iz)*  A  condition  sufficient  to  insure  the  accuracy  of 
(2.2)  is 

|  ^2  *2  |  >>  ^1  tl  =  tQ  (2.4) 

where  e2  is  given  by 

42  =  t*  -  j  o2/u  ,  e*  >  0  (2.5) 

and  and  a^  are  the  permittivity  and  permeability  of  free 

space.  In  (2.5),  o2  is  the  conductivity  ofthe  half-space  M2  and  a 
harmonic  time  variation  of  the  form  exp(jwt)  has  been  assumed. 

When  (2.4)  is  satisfied,  the  fields  in  M2  are  constrained 
to  propagate  along  the  z  axis  (in  the  -z  direction)  like  plane 
waves,  and  the  validity  of  (2.2)  is  assured.  The  components  E2x> 
H2  ,  E2  ,  and  H2  satisfy  the  one-dimensional,  homogeneous  wave 

y  y  x 

equation  with  solutions  of  the  form 

E 2 x ( x , y , z )  =  E2x(x,y,0)  exp(i2  z) 

H2y(x,y,z)  =  H2y(x,y,0)  exp(i2  z)  (2.6) 

E2  (  x  ,  y ,  z  )  =  E2  (  x  ,  y  ,  0  )  exp(-»2  z) 

H2x(x,y,z)  =  H2x(x,y,0)  exp(i2  z) 

where  E2  (x,y,0),  H2  (x,y,0),  E2  (x,y,0),  and  H2  (x,y,0)  are  the 
x  y  y  x 

the  tangential  components  of  the  field  transmitted  into  M2  at  the 
interface  z  =  0  due  to  the  fields  in  Ml  and  the  propagation  con¬ 
stant  i2  is  given  by 

f2  =  jw(p2  1 2 ) 1 /2  .  (2,7) 

When  E2  ( x , y , 0 ) ,  H2  (x,y,0),  E2  (x,y,0),  and  H2  (x,y,0)  are 
x  y  y  x 


known,  equation  (2.6)  can  be  used  to  compute  the  fields  at  any 
point  in  M2. 

Equation  (2.2)  is  frequently  referred  to  as  the  Leontovich 
5 

Boundary  Condition  although  it  was  apparently  used  before  Leon- 

6  7  8 

tovich  by  Rytov,  Alpert,  and  Feinberg0  during  World  War  II  in 
their  work  on  ground  wave  propagation.  Since  then,  these  and 
other  investigators  have  shown  that  (2.2)  is  not  limited  in  ap¬ 
plication  to  half-spaces  but  can  be  applied  directly  to  more 
complicated  structures.  For  example,  if  the  half-space  M2  is  re¬ 
placed  by  a  sheet  of  the  same  material  with  a  uniform  thickness  d 
as  shown  in  figure  3,  then  the  fields  inside  the  sheet  still  sat¬ 
isfy  the  one-dimensional  wave  equation  and  propagate  like  plane 
waves  parallel  to  the  z  axis  provided  (2.4)  remains  valid.  The 
IBC  is  again  applicable  at  z  =  0  in  the  form  given  by  (2.2)  if  an 
additional  condition  is  satisfied: 

6g  <  d  (2.8) 

where 

6g  =  1/u  |  Im( (p2e2) 1/2)  |  (2.9) 

is  the  skin  depth  of  the  sheet  (M2.)  The  skin  depth  is  a  measure 

5 

M.  A.  Leontovich,  Investigation  of  Radiowave  Propagation,  Part 
II,  Moscow:  Academy  of  Sciences  (1948). 

6  S.  M.  Rytov,  J.  Exp.  Theor.  Phys .  USSR,  10  (1940), 180. 

7  I.  L.  Alpert,  J.  Tech.  Phys.  USSR,  10  (1940),  1358. 

8  E.  L.  Feinberg,  J.  Phys.  USSR,  8  (1944),  317. 
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Figure  3.  A  planar  sheet  of  uniform  thickness  d 
separating  media  where  an  impedance  boundary  condition 
is  satisfied  at  z  =  0  by  virtue  of  the  relations 
1  tl  <<  |^2  t2|  and  6g  <  d.  (Reproduced  from  ref  1.) 
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of  the  rate  at  which  the  magnitude  of  the  field  decreases  as  it 
propagates  in  M2.  Since  fields  decrease  in  amplitude  by  a  factor 


of  e 


0.37  =  8.5  dB  while  traveling  a  distance  equal  to  one 


skin  depth,  condition  (2.8)  means  that  fields  making  a  round  trip 

from  z=0toz=-d  and  back  to  z  =  0  will  be  reduced  by  at 

_2 

least  a  factor  of  e  =  0.14  =  17  dB .  This  condition  is  necessary 
to  prevent  fields  reflected  at  z  =  -d  from  reaching  z  =  0  in 
sufficient  strength  to  interfere  with  El  and  Hi  at  the  interface 
and  render  (2.2)  inaccurate.  This  also  means  that  internally  re¬ 
flected  fields  at  z  =  0  can  be  neglected  and  that  the  fields  at  z 
=  0  traveling  into  M2  are  again  given  by  (2.6).  In  other  words, 

(2.8)  eliminates  standing  waves  in  the  sheet,  and  a.lows  one  to 
represent  the  fields  propagating  into  the  sheet  at  z  =  0  as  sim 
pie  traveling  waves  just  as  in  the  case  of  the  infinite  half¬ 
space.  Condition  (2.8)  is  easily  satisfied  by  most  electromag¬ 
netic  shields  composed  of  flat  sheets  since  these  structures  are 
designed  to  reduce  the  fields  reaching  z  =  -d  by  far  more  than 
8.5  dB.  In  fact,  it  will  usually  be  found  that  quality  shields 

satisfy  the  much  more  stringent  condition  6  <<  d  at  most  fre- 

s 

quencies  of  interest. 

If  a  rectangular  slot  is  cut  through  the  uniform,  homogen¬ 
eous  sheet  in  figure  3  and  the  rectangular  volume  of  the  slot  is 
filled  by  a  material  with  permeability  2'  and  permittivity  t2’, 
then  the  sheet  is  rendered  locally  inhomogeneous  and  anisotropic 
-  inhomogeneous  because  the  effective  permeability,  permittivity, 


and  impedance  of  the  slot  differ  from  those  of  the  solid  sheet 
and  anisotropic  because  these  quantities  depend  on  the  orienta¬ 
tion  of  the  slot.  If  the  slot  is  oriented  as  shown  in  figure  (a), 

then  El  and  HI  will  see  a  different  impedance  looking  into  the 
x  y 

slot  than  will  El  and  HI  .  With  these  two  impedances  denoted  n 2 

y  x  x 


and  q2,  a  generalized  IBC  can  be  written 

El  =  -  q2  HI  ,  El  =  i?2  HI  ,  at  z  =  0 
x  x  y  y  y  x 


(2.10) 


where  q 2  and  q 2  are  functions  of  x  and  y,  p 2'  and  t2'  ,  and  the 
x  J 

dimensions  of  the  slot  are  a  and  p.  A  similar  IBC  can  be  written 
for  the  vertical  slot  shown  in  figure  4. 

A  condition  sufficient  to  insure  that  the  fields  in  the 
slot  E2^,  H2^.,  E2^,  and  H2^  propagate  parallel  to  the  z  axis  in 
the  manner  of  plane  waves  is 


<  <  Z 


<  W 


(2.11) 


Since  q 2  and  q2  can  be  defined  in  terms  of  effective  slot 
x  y 

permeabilities  and  permittivities, 


<72  =  <„2/e2 


<72  =  <*$/*$) 


(2.12) 


where  p2  ,  p2  >  Condition  (2.10)  then  implies 


I  ‘i  I 


|  y2  «2  I 

I  ey  y  | 


>  >  u  t 
o  o 


(2.13) 


which  is  completely  analogous  to  (2.4).  The  fields  in  the  slots 


satisfy  separate  homogeneous  wave  equations  corresponding  to  q2 
and  q 2.  The  solutions  to  these  equations  give  the  following 

y 

expressions  for  the  slot  fields  at  z  =  0  propagating  in  the  -z 
direction : 

E2^(x,y,z)  =  E2^(x,y,0)  exp(i2x  z) 


H2 ' (x ,  y ,  z  )  =  H2  ‘  ( x , y , 0 )  exp(i2  z) 
y  y  A 


(2.14) 


E2  '  ( x , y , z )  =  E2 ' ( x , y , 0 )  exp(i2  z) 
y  y  y 

H2^(x,y,z)  =  H  2  ^  ( x ,  y ,  0  )  exp(i2y  z) 
where  the  propagation  constants  are  given  by 
ir2  =  jw  p2  /  q2 

72  =  jw  p2  /  ,2  .  (2-15) 

y  *y  'y 

For  air-filled  slots  with  dimensions  that  are  small  com¬ 
pared  to  the  wavelength  of  El, HI,  the  slot  permeabilities  p2  and 

p2.  are  real  and  equal  to  the  permeability  of  free  space  u  and 
J  O 

the  slot  impedances  are 


nl  *  Lx 


”1  *  Ly 


(2.16) 


where  L  ,  L  >0  are  inductances.  In  this  case,  (2.15)  reduces 
x  y 

to  72  =  /L^  and  72  =  /Ly  .  Equation  (2.14)  becomes 

E2^(x,y,z)  =  E2^(x,y,0)  exp(z 


H2 ' ( x , y , z )  =  H2 ' ( x , y , 0 )  exp(z  p  /L  ) 
y  y  u  a 

E2^(x,y,z)  =  E 2 ^ ( x , y , 0 )  exp(z  PQ/Ly) 


(2.17) 


H2^(x,y,z)  =  H2^(x,y,0)  exp(z  P0/Ly! 


The  latter  show  that  the  slot  acts  like  a  waveguide  below  cutoff 
resulting  in  an  exponential  decay  of  the  slot  fields  for  z  <  0. 
The  rates  of  decay  depend  primarily  on  the  dimensions  of  the  slot 


a  and  p .  As  a  and  p  decrease,  L  and  L  also  decrease,  and  the 

x  y 

rate  of  decay  increases.  The  air-filled  slot  is  a  structure  of 

g 

considerable  importance  in  shielding  theory  since  Jarva  has 
shown  that  it  can  be  used  as  a  working  model  for  the  most  common 
types  of  defects  (discontinuities)  that  occur  in  practical 
shields 

The  slotted  sheet  like  the  homogeneous  sheet  must  satisfy 

an  additional  condition  involving  skin  depth  and  sheet  thickness 

in  order  to  insure  that  (2.10)  is  valid  at  the  interface  between 

Ml  and  the  slot.  Here  there  are  two  skin  depths,  5  and  6  , 

x  y 

which  in  the  case  of  the  air-filled  slot  can  be  written  as 
follows : 


6  =  L  /  u 

x  x  o 


0  =  L  /  u 

y  y  o 


(2.19) 


and  the  condition  (analogous  to  (2.8)  is 


<  d 


(2.20) 


This  condition  insures  that  the  impedance  at  z  =  0  is  unaffected 
by  reflections  at  z  -  -d.  When  (2.13)  and  (2.20)  are  satisfied, 
the  IBC  (2.10)  is  a  valid  local  boundary  condition  at  the 
interface  between  Ml  and  the  slot,  and  the  fields  in  the  slot 


W.  Jarva,  IEEE  Trans.  EMC,  EMC-12  (1970),  12, 
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traveling  away  from  the  interface  are  plane  waves  given  by 
(2.14).  If  a  sheet  has  more  than  one  slot,  then  (2.10)  can  be 
applied  at  each  slot  provided  the  appropriate  conditions  are 
satisfied  and  the  impedances  of  the  slots  can  be  determined. 

2.2  Shielding  by  Continuous  Structures  Composed  of  Good 
Conductors 

IBC’s  were  developed  originally  to  simplify  external 
scattering  and  propagation  problems  in  which  the  sole  objective 
is  to  compute  the  fields  in  Ml: 

eT  =  eTs  +  Elr 

(2.21 ) 

Hi  =  Hi®  +  hT1' 

where  El  and  HI  are  the  fields  generated  by  the  source  and  El 1 

and  Hlrare  the  reflected  and  diffracted  fields  in  Ml  due  to  M2. 
In  most  cases,  the  source  fields  are  known  functions  of  position 
in  Ml  and  the  external  problem  reduces  to  that  of  determining  the 
reflected  and  diffracted  fields.  An  IBC  applied  at  the  interface 
between  Ml  and  M2  simplifies  such  a  problem  by  decoupling  the 
fields  in  Ml  from  those  in  M2.  This  allows  one  to  solve  the  ex¬ 
ternal  problem  without  the  necessity  of  solving  the  internal 

problem  for  E2  and  H2 .  However,  in  many  applications,  such  as 
underground  communication  and  electromagnetic  shielding,  the 
internal  problem  is  of  equal  or  greater  importance  than  the 
external  problem,  and  it  is  natural  to  ask  if  an  IBC  can  be  used 


to  simplify  the  problem  of  computing  E2  and  H2.  It  was  shown  by 
Monroe/  that  an  IBC  can  indeed  be  used  to  obtain  appoximate 
solutions  to  certain  internal  problems  by  relatively  simple  means 
and  that  the  latter  lead  naturally  to  formal  solutions  to  a  gen¬ 
eral  class  of  shielding  problems. 

Monroe  considers  the  problem  of  computing  fields  inside  the 

generalized  structure  shown  in  figure  5  when  it  is  illuminated  by 

an  arbitrary  source  S.  The  structure  consists  of  flat  sheets  ESj  , 

ES_,  . . . , ES  of  uniform  composition  and  thickness  d  attached  con- 
c  m 

tinuously  along  their  edges.  This  structure  includes  both  open- 
and  closed-form  shields  as  special  cases  where  M2  corresponds  to 
the  shield  itself  and  M3  is  the  shielded  volume.  By  specifying 
the  number,  shape,  and  dimensions  of  the  sheets,  one  can  generate 
the  half-space  and  the  infinite  sheet  described  earlier  and  the 
rectangular  parallelepiped  among  many  other  useful  arrangements 
of  flat  surfaces. 

The  composition  of  the  sheets  is  assumed  to  fall  into  the 
class  of  materials  referred  to  as  good  conductors,  which  is  char¬ 
acterized  by  the  relation: 

w  t*  <<  o2  (2.22) 

where  t*  is  the  real  part  of  the  complex  permittivity  (2.5)  and 
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(a)  Sideview 


z  ■< 


(b) 


Front  View 


a 2  is  the  conductivity.  For  these  materials,  r}2  and  t 2  reduce  to 
r}2  =  (jw  p2/o2)1/2  (2.23) 

i2  =  (jw  ij2o2  ) 1  /2  .  (2.24) 

Since  all  primary  materials  used  in  shield  contruction  satisfy 
(2.22)  at  all  frequencies  of  interest,  and  since  (2.22)  is  con¬ 
sistent  with  (2.4),  it  is  clear  that  Monroe  did  not  impose  a  se¬ 
rious  limitation  on  the  theory  by  restricting  it  to  structures 
composed  of  good  conductors.  On  the  other  hand,  by  accepting  this 
limitation,  one  is  able  to  simplify  the  theory  using  the  well- 
known  fact^  that  the  net  tangential  component  of  the  magnetic 
field  at  the  surface  of  a  good  conductor  is  approximately  equal 
to  twice  the  corresponding  component  of  the  incident  (source) 
field. 

The  objective  is  to  obtain  formal  expressions  for  the  elec¬ 
tric  and  magnetic  fields  transmitted  into  the  shielded  volume  M3 
at  any  point  P'(x,y,-d)  on  the  inside  surface  of  ES^  due  to  field 
from  S  incident  on  the  outside  surface  of  ESj  (z  =  0)  where  the 
fields  are  referred  to  a  rectangular/cylindrical  coordinate  sys¬ 
tem  that  is  consistent  with  the  system  used  to  define  the  IBC. 
The  origin  of  this  system  is  located  on  the  outside  surface  of 
ESj  and  the  positive  z  axis  points  away  from  M3.  For  convenience, 


R.  B.  Adler,  L.  J.  Chen,  and  R.  M.  Fano,  Electromagnetic 
Energy  Transmission  and  Radiation,  John  Wiley  and  Sons,  Inc. 
N.Y.  (1960)  p  432. 
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the  source  is  located  on  the  z  axis,  the  origin  is  located  at  the 
geometrical  center  of  ES^,  and  the  x  and  y  axes  are  oriented 
parallel  to  the  edges  of  ES^  as  shown  in  figure  5. 

The  transmitted  fields  at  P’(x,y,-d)  are  computed  in  terms 
of  the  tangential  components  of  the  source  magnetic  field  inci¬ 
dent  on  ESj  by  a  five-step  process:  In  the  first  step,  the  IBC  is 
used  together  with  standard  boundary  conditions  to  show  that  the 
x  and  y  components  of  the  electric  and  magnetic  fields  at  z  =  0 
in  ESj  can  be  written  in  terms  of  the  corresponding  components  of 
the  magnetic  field  in  Ml.  That  is, 

E2  (x,y,0)  =  -  q2  HI  (x,y,0) 

x  y 


H2  ( x , y , 0 )  =  HI  ( x , y , 0 ) 

y  y  (2.25) 

E2  ( x , y , 0 )  =  q2  HI  (x,y,0) 
y  ^ 

H2x(x,y,0)  =  HI x ( x , y , 0 )  . 

The  second  step  uses  the  previously  noted  property  of  good  con¬ 
ductors  to  write  HI  (x,y,0)  and  HI  (x,y,0)  in  terms  of  the  cor- 

x  y 

responding  components  of  the  source  magnetic  field: 


Hlx(x,y,0)  s  2  HI “ ( x , y , 0 ) 


Hly( x , y , 0 )  s  2  Hly(x,y,0) 


(2.26) 


and  subsitutes  these  expressions  in  (2.25)  to  obtain 


E2  ( x , y , 0 )  =  -  2  rj2  Hl°(x,y,0) 
x  y 

H2  ( x , y , 0 )  =  2  Hl®(x,y,0) 

y  <7 

E2  (x,y,0)  =  2  17 2  Hls(x,y,0) 


(2.27) 


CT.V.  AV 


H2x(x,y,0)  =  2  Hl"(x,y,0). 

The  next  step  is  to  substitute  (2.27)  into  (2.6)  to  obtain  ex¬ 
pressions  for  the  plane  wave  fields  in  M2  traveling  toward  M3  in 
terms  of  the  tangential  components  of  the  source  magnetic  field 
incident  on  ES^: 

E2  ( x ,  y ,  z )  =  -  2  q2  Hl®(xfy,0)  exp(i2  z) 
x  y 


H2  ( x , y , z )  =  2  HI  ( x , y , 0 )  exp(i2  z) 

(2.28) 

E2  ( x , y , z )  =  2  q2  Hl®(x,y,0)  exp(i2  z) 

y  x 

H2x(x,y,z)  =  2  Hlx(x,y,0)  exp(i2  z) 

where  q2  and  i2  are  given  by  (2.23)  and  (2.24).  Next,  the  plane 
wave  fields  incident  normally  on  the  interface  between  M2  and  M3 
are  obtained  by  evaluating  (2.28)  at  z  =  -  d: 

E2  ( x , y , -d )  =  -  2  q2  Hl®(x,y,0)  exp(-  i2  d) 
x  y 

H2  ( x , y , -d )  =  2  Hl®(x,y,0)  exp(-  i2  d) 

(2.29) 

E2  ( x  ,  y ,  -d )  =  2  rj2  Hl®(x,y,0)  exp(-  t2  d) 

y  ^ 

H2x ( x , y , -d )  =  2  q2  Hlx(x,y,0)  exp(-  i2  d). 

These  generate  reflected  fields  in  M2  traveling  in  the  +z  direc¬ 
tion  and  transmitted  fields  in  M3  traveling  in  the  -z  direction. 
The  final  step  is  to  compute  the  transmitted  field  at  z  =  -  d 
from  (2.29).  This  is  done  by  multiplying  (2.29)  by  appropriate 
plane  wave  trnsnsmiss ion  coefficients.  The  result  is: 

E3x(x,y,-d)  =  -  2  772  Tg  Hl®(x,y,0)  exp(-  i2  d) 


(2.29) 


H3  (x.y.-d)  =  2  T  Hl"(x,y,0)  exp(-  ?2  d) 

j  n  jr 

E3  (x,y,-d)  =  2  q2  T_  Hl®(x,y,0)  exp(-  i2  d) 

y  £  x 


(2.30) 


J  —  - 


H3x(x,y,-d)  =  2  Th  Hl®(x,y,0)  exp<-  i2  d) 

where  T„  and  Tu  are  the  transmission  coefficients  for  electric 
.h  n 

and  magnetic  plane  wave  fields  incident  normally  on  a  planar  sur¬ 
face  separating  two  dissimilar  media.  These  coefficients  are 
2 

given  by: 

T£  =  2  q3/(q2  +  q3) 

(2.31) 

TH  =  2  q2/(q2  +  q3) 

where  q2  is  as  previously  defined  and  q3  is  the  impedance  at 
z  =  -d  looking  into  M3. 

The  preceding  are  formal  expressions  for  the  principal  com¬ 
ponents  of  the  electric  and  magnetic  fields  at  the  surface  of  the 
shielded  volume  (z  =  -d )  for  the  generalized  structure  in  figure 
5.  To  obtain  explict  expressions  for  a  particular  shielding  prob- 

g 

lem,  it  is  necessary  to  determine  the  source  fields  Hlx(x,y,0), 
Hl®(x,y,0),  and  the  impedance  q 3.  For  example,  if  S  is  an  antenna 

g 

located  a  specified  distance  from  ESj,  then  HI  (x,y,0)  and 

Hl®(x,y,0)  can  be  computed  using  standard  methods  from  antenna 

theory.  This  was  done  in  Monroe^  for  Hertzian  dipoles  and  small 

S  8 

rectangular  loop  antennas.  However,  if  HI  (x,y,0)  and  HI  (x,y,0) 

x  y 

2 

S.  A.  Schelkunoff,  Electromagnetic  Waves,  Van  Nostrand,  Prince¬ 
ton,  NJ  (1943). 
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are  known  functions  or  can  be  estimated  from  measurements,  then 
these  quantities  can  be  used  directly  in  (2.30),  and  nothing  more 
need  be  specified  about  the  source.  This  will  be  useful  in  many 
cases  where  the  source  (or  sources)  generating  Hl®(x,y,))  and 

g 

HI  (x,y,0)  is  partially  or  completely  unknown.  Examples  of  the 

y 

latter  include  lightning  strokes  and  NEMP  (nuclear  electromagnet¬ 
ic  pulse).  By  the  same  token,  it  follows  that  two  or  more  sources 
generating  the  same  magnetic  field  components  tangent  to  ES^  will 
produce  the  same  electric  and  magnetic  fields  in  M3. 

To  determine  q 3,  it  is  necessary  to  specify  the  geometry  of 
the  structure.  For  open  structures,  q 3  can  usually  be  approxi¬ 
mated  by  the  wave  impedance  of  the  source.  This  was  done  by  Mon¬ 
roe  for  the  case  of  an  infinite  flat  sheet  with  free  space  on 
either  side  exposed  to  Hertzian  dipoles.  In  this  case,  the  ap¬ 
proximation  is  justified  by  the  evident  fact  that  the  structure 
of  the  field  must  be  nearly  the  same  on  both  surfaces  of  the 
sheet  since  Ml  and  M3  are  both  of  infinite  extent.  However,  the 
same  is  not  true,  in  general,  if  M2  is  closed.  In  this  case,  the 
structure  of  the  fields  in  M3  usually  has  no  simple  relationship 
to  that  of  the  source  field  incident  on  ES^.  Instead,  the  struc¬ 
ture  of  the  fields  in  M3  is  determined  by  the  geometry  of  the  en¬ 
closure  and  by  its  size  relative  to  the  free  space  wavelength  of 
the  source  field.  If  the  cross  section  of  M3  transverse  to  the  z 
axis  is  uniform  (does  not  change)  from  z  =  -  d  to  z  =  -L  where  L 


denotes  the  surface  of  the  sheet  opposite  ES^  ,  then  M3  is  equiv¬ 
alent  to  the  interior  of  a  section  of  waveguide  closed  at  both 
ends.  In  this  case,  q 3  will  be  determined  by  one  or  more  of  the 
doubly  infinite  set  of  waveguide  modes  that  this  structure  can 
support.  In  Monroe,*  the  principal  features  of  fields  transmitted 
to  the  interior  of  a  rectangular  enclosure  were  approximated  at 
the  point  z  =  -d  using  a  single  rectangular  waveguide  mode-  the 
dominant  TE^q  mo<ie  .  In  section  3,  this  approach  is  extended  by 
expanding  the  fields  as  doubly  infinite  series  over  all 
waveguide  modes  and  applying  equations  (2.29)  and  (2.30)  to  each 
mode  at  z  =  0 . 

2 . 3  Shielding  by  Discontinuous  (Slotted)  Stuctures 
Composed  of  Good  Conductors. 

If  the  structure  in  figure  5  is  modified  by  cutting  n 

narrow  rectangular  slots  through  ES ,  at  various  locations  (x^y1) 

1  c '  c 

for  i  =  l,2,3,....n  ,  where  x1  and  y1  are  the  x  and  y  coordinates 

c  c 

of  the  center  of  the  i-th  slot,  then  the  continuous  structure  is 
transformed  into  one  with  multiple  discontinuities.  Using  (2.10), 
and  (2.14)  and  a  5  step  process  analogous  to  the  one  used  for  the 
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continuous  structure,  Monroe  obtained  expressions  for  the  fields 
at  the  center  of  the  i-th  slot  incident  on  M3  at  z  =  -  d  and  the 
fields  transmitted  to  M3  at  the  same  point.  These  expressions 
were  then  used  to  describe  the  principal  features  of  the  fields 
at  the  center  of  a  slot  in  one  wall  of  a  rectangular  enclosure.  A 
similar  process  yields  the  following  expressions  for  the  fields 
incident  on  the  surface  of  the  i-th  slot  at  z  =  -d. 

E2*(x,y,-d)  =  -  2  rjb  Hl®(x,y ,0)  fe*(x,y,-d) 

x  exp ( -jwd  v\/n\) 


H2*(x,y,-d)  =  2  Hl®(x,y,0)  fh*(x,y,-d)  exp(-jwd 

y  y  y  XX 


(2.32) 


E2*(x,y,-d)  =  2  Hl®(x,y,0)  fe*(x,y,-d) 

y  y  x  y 


x  exp(  -  jwd  fj$/np 

H2 * ( x , y , -d )  =  2  Hl®(x,y,0)  fh*(x,y,-d)  exp(-jud  v\/n\) 
x  xx  y  y 

where  all  quantities  are  as  previously  defined  except  the  form 

factors  fe1,  fh1,  fe1.  and  fh1  which  describe  the  spatial  varia- 
x  y  y  x 

tion  of  the  incident  electric  and  magnetic  field  components  at  z 
=  -  d  due  to  the  i-th  slot.  The  fields  transmitted  to  M3  from  the 
i-th  slot  at  z  =  -d  are  obtained  from  (2.32)  by  applying  appro¬ 
priate  transmission  coefficients.  That  is, 

E3*(x,y,-d)  =  T£*  E2*(x,y,-d) 


H3*(x,y,-d)  =  TH*  H2*(x,y,-d) 

*  J  ¥ 

E3*(x,y,-d)  =  TE*  E2*(x,y,-d) 
H3*(x,y,-d)  =  TH*  H2*(x,y,-d) 


(2.33) 
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ARBITRARY,  EXTERNAL,  TIME  HARMONIC,  ELECTROMAGNETIC  FIELDS 


The  formal  expressions  in  section  2  can  be  adapted  to  give 
the  fields  that  penetrate  a  rectangular  enclosure  by  a  process 
analogous  to  the  classical  plane  wave  expansion  technique  where 
an  electromagnetic  field  is  represented  as  the  sum  of  an  infinite 
series  of  plane  waves**.  Since  any  source  field  can  be  repre¬ 
sented  in  this  way,  plane  wave  expansions  have  often  been  used  to 
simplify  external  problems  (figure  1)  by  reducing  the  fields  in 
Ml  to  the  sum  of  the  source  field  and  a  series  of  plane  wave 
reflections  from  M2.  Indeed,  the  same  process  could  be  applied  to 
the  shielding  problem  by  computing  the  plane  waves  transmitted 
from  Ml  to  M2.  However,  the  impedance  boundary  condition  makes 
this  unnecessary  since  it  gives  the  plane  wave  fields  in  M2 
directly  in  terms  of  the  source  fields  tangent  to  the  interface 
between  Ml  and  M2. 


A  different  expansion  is  needed  for  the  shielding  problem: 
one  that  can  be  used  to  relate  the  plane  wave  fields  incident  at 
the  M2/M3  interface  (equations  (2.29)  and  (2.32))  to  the  trans¬ 
mitted  fields  at  that  interface  and  to  the  fields  in  the  interior 
of  M3.  Since  a  rectangular  enclosure  is  equivalent  to  a  section 


**P.  C.  Clemmow,  The  Plane  Wave  Representation  of  Electromagnetic 
Fields,  Pergamon  Press,  Oxford  (1966). 
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of  rectangular  waveguide,  the  natural  choice  here  is  an  expansion 
in  terms  of  the  mode  functions  for  this  structure.  Expanding 
(2.29)  and  (2.32)  into  series  composed  of  mode  functions  express¬ 
es  each  of  these  plane  wave  fields  as  the  sum  of  an  infinite  ser¬ 
ies  of  plane  waves  each  of  which  corresponds  to  a  single  wave¬ 
guide  mode.  In  effect,  the  expansions  replace  single  plane  wave 
fields  by  an  infinite  number  of  elementary  fields  that  are  both 
plane  waves  and  mode  functions.  Each  of  these  elementary  fields 
is  incident  normally  on  the  M2/M3  interface  and  each  generates  a 
reflected  field  in  M2  traveling  in  the  +z  direction  and  a  trans¬ 
mitted  field  in  M3  traveling  in  the  -z  direction.  The  transmitted 
field  at  the  interface  can  be  computed  in  terms  of  the  incident 
elementary  field  by  using  an  appropriate  plane  wave  transmission 
coefficient  in  the  same  manner  that  equations  (2.30)  and  (2.33) 
were  obtained  from  (2.29)  and  (2.32).  Since  the  transmitted  field 
at  the  interface  is  also  a  wavequide  mode  for  the  enclosure,  it 
can  be  used  to  obtain  the  field  at  any  point  inside  M3  due  to  the 
incident  elementary  field  at  the  interface  by  representing  the 
field  in  M3  by  the  same  mode  and  applying  boundary  conditions  at 
the  walls  and  interface.  The  advantage  to  this  approach  is  that 
the  total  field  at  any  point  in  M3  can  be  obtained  by  summing  the 
contributions  at  that  point  due  to  all  the  elementary  fields  in 
the  original  expansion.  This  sum  is  a  complete  solution  to  the 
shielding  problem  since  it  is  the  field  inside  M3  due  to  the 


external  source. 


3.1  Normal  Mode  Functions  for  an  Enlosure  Composed  of  a 
Good  Conductor 

The  shell  structure  shown  in  figure  6  is  the  ideal  ver¬ 
sion  of  the  most  common  practical  electromagnetic  shield.  It  is  a 
special  case  of  the  generalized  structure  (figure  5)  in  the  form 
of  a  rectangular  parallelepiped  with  walls  composed  of  a  good 
conductor  as  defined  in  section  2.  The  walls  enclose  a  volume  of 
free  space  (M3)  with  arbitrary  dimensions  A,  B,  and  C.  As  in 
figure  5,  we  consider  fields  from  an  electromagnetic  source  S 
entering  M3  through  the  front  wall  of  the  enclosure,  that  is, 
through  the  wall  directly  exposed  to  the  source.  (If  two  or  more 

of  the  walls  are  exposed  to  the  same  source,  or  to  a  different 
source,  then  the  techniques  developed  here  can  be  applied  to  each 
wall  in  turn,  and  the  total  field  at  any  point  in  M3  can  be  ob¬ 
tained  by  adding  contributions  from  all  the  walls.)  The  rectang¬ 
ular/cylindrical  coordinate  system  with  origin  at  the  geometric 
center  of  the  outside  surface  of  the  front  wall  is  oriented  with 
x  and  y  axes  parallel  to  the  edges  of  the  wall.  With  the  +z  axis 
pointing  outward,  the  inside  surface  of  the  front  wall  lies  in 
the  z  =  -d  plane  where  d  is  the  wall  thickness,  the  inside  sur¬ 
face  of  the  back  wall  lies  in  the  z  =  -(d  +  C)  plane,  and  the 

outside  surface  of  the  back  wall  lies  in  the  z  =  -(2d  +  C)  plane. 

Similarly,  the  inside  surfaces  of  the  side  walls  lie  in  the  x  =  - 

A/2  and  y  =  -B/2  planes. 
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(a)  Side  View 


( 


4 


(b)  Front  View 

Figure  6. 
continuous 
parallelep 
dimensions 


With  this  arrangement  M3  is  defined  as  follows 

-  A/2  <  x  <  A/2 

M3:  -  B/2  <  y  <  B/2  (3.1) 

-  (d+C)  <  z  <  -d 

and  it  is  seen  that  for  fields  traveling  into  the  enclosure  along 

paths  parallel  to  the  z  axis,  the  M2/M3  region  is  a  rectangular 

waveguide  driven  at  z  =  -  d  and  loaded  at  z  =  -  (d  +  C)  by  a 

sheet  with  finite  conductivity 

The  fields  in  M2/M3  can  be  represented  most  conveniently  in 

12 

terms  of  the  normalized  TE  and  TM  mode  functions  for  the  equiv¬ 
alent  waveguide.  These  functions  are  solutions  to  Maxwell’s  equa¬ 
tions  satisfying  appropriate  boundary  conditions  at  the  side 
walls  x  =  +  A/2  and  y  =  +  B/2.  The  two  types  of  modes  are  dis¬ 
tinguished  by  the  fact  that  TE  modes  have  no  electric  field  com¬ 
ponent  in  the  z  direction  and  TM  modes  have  no  magnetic  component 
in  the  z  direction.  In  vector  form,  the  TE  mode  functions  refer¬ 
red  to  the  coordinate  system  in  figure  6  are 


er.m(x»y) 

nm 


=  K 


nm 


nm .  . 

e  (  x  ,  y  )  l 
x  x 


nm ,  . 

e  (  x  ,  y )  l 

y  y 


where  i  and  i  are  unit  vectors, 
x  y 

e™(x,y)  =  (mn/B)  cos(nn(x  +  A/2)/A)  sin[mn(y  + 


ey  (x,y)  =  -  (nn/A)  sin[nn(x  +  A/2)/A)  cos[nm(y 
1  2 

R.  E.  Collin,  Field  Theory  of  Guided  Waves, 
York  (1960). 


(3.2) 


B/2)/B] 


+  B/2)/B] 

McGraw-Hill,  New 
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K  =  ( nn/A )  +  (mn/B)‘  [e  t  /AB]1"6 

nm  on  on 


e  and  t 
on  om 


on 


1 

2 


are  Neuman’s 

when  n  =  0 
when  n  >  0 


f  actor : 


and  all  other  quantities  are  as  defined  previously.  These  func¬ 
tions  give  the  spatial  structure  of  the  electric  fields  associ¬ 
ated  with  the  TE  modes.  The  structure  of  the  transverse  compon¬ 
ents  of  the  magnetic  fields  associated  with  the  TE  modes  is  given 
by 


h„_,  ( x  ,  y )  =  K _ reym(  x  ,  y  )  i^  +  e^m(x,y)  i 


nm 


nm 


(3.4) 


where  all  quantities  are  as  defined  previously.  Similar  expres¬ 
sions  for  TM  modes  can  be  obtained  by  simple  transformations  of 
(3.2)  and  (3.4). 


When  e  (x,y)  and  h  (x,y)  are  multiplied  by  factors  of  the 
nm  nm  r 


form  exp(+  rnmz)  where  fnm  is  the  propagation  constant  for  the 


TE  mode,  the  products  are  also  solutions  to  Maxwell's  equa- 
nm 


tions.  That  is,  equations  (3.5)  and  (3.6) 


ft  —  nm 

e  ( x , y , z )  =  e  (x,y)  exp(+r  z) 


nm 


nm 


(3.5) 


h*m(x,y,z)  =  hnm(x,y)  exp(+  rnmz) 


where 

.nm 


1/2 


:c  ) 


-l 


(3.6) 


r*.' 


(3.7) 


p.v.v.v 


»7»’."/  «/  «,* 


and  w 


nm 


c 

nm 


is  the  cutoff  frequency 

1/2 


(n/ 


«c  =  nCo|(n/A)2  +  (m/B)21 


are  solutions  to  Maxwell’s  equations  at  all  points  in  the  wave¬ 
guide.  When  the  exponential  carries  a  "+"  sign,  (3.5)  represents 

TE  mode  fields  traveling  in  the  -z  direction:  conversely,  a 
nm 

sign  represents  fields  traveling  in  the  +z  direction.  Furthermore 


since 

Maxwell’s  equations 

are 

linear,  fields 

in 

E31 

i  in  i  v  rrnin  <  v 

( x , y , z )  =  E3  ( x , y ) 

rnm 

fe 

(  z ) 

H37 

( x  ,  y ,  z  )  =  H3nm  ( x  ,  y  ) 

„nm 

fh 

(  z )  , 

where 

E31 

un  .  .  .nm  nm ,  . 

(x,y)  =  $x  ex  ( x , y ) 

ix 

.nm  nm ,  . 

+  (y  ey  ( X , y ) 

* 

i 

y 

H37 

im,  .  nm  nm.  . 

(  x  ,  y  )  =  {  e  (  x  ,  y  ) 
x  y 

4x 

.nm  nm,  . 

+  Cy  ex  ( x , y ) 

a 

i 

y 

(3.8) 


(3.9) 


F”m(z)  =  Cl  exp(+rnmz)  +  C2  exp(-rnmz) 

F„m(z)  =  C3  exp(+/"nmz)  +  C4  exp(-rnmz) 
H 


and  $nm,  ?nm,  {nm,  {nm,  Cl,  C2,  C3,  and  C4  are  constants,  are 
x  y  x  y 

also  solutions  to  Maxwell’s  equations.  These  expressions  repre¬ 
sent  waveguide  fields  as  sums  of  fields  traveling  in  the  +z  and 
-z  directions,  such  as  would  occur  if  fields  in  the  waveguide 
were  reflected  between  a  driver  and  a  termination.  By  applying 
appropriate  boundary  conditions  at  the  locations  of  the  driver 
and  the  termination,  say  z  =  -d  and  z  =  -(d  +  C),  the  constants 


in  (3.9)  can  be  computed.  With 


,nra 


(z)  and  F^m(z)  determined, 
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(3.8)  becomes  the  correct  form  for  the  mode  fields  in  a  finite 
section  of  waveguide  driven  at  z  =  -d  and  terminated  at  z  =  -d 
C.  That  is,  (3.8)  becomes  the  correct  form  for  the  mode  fields  in 


Although  a  determination  of  F”m(z)  and  Fi]m(z)  is  ultimately 

is  H 

based  on  applying  boundary  conditions  at  the  specfied  locations, 
we  can  obtain  these  factors  directly  without  the  lengthy  calcu¬ 
lations  required  in  such  a  procedure  by  employing  the  fact  that 
the  mode  functions  (3.8)  satisfy  standard  transmission  line  equa¬ 
tions.  This  fact  establishes  an  exact  mathematical  analogy  be¬ 
tween  the  mode  functions  and  solutions  to  transmission  line  prob¬ 
lems  that  allows  us  to  transform  the  latter  into  the  former  using 
simple  substitutions  such  as 

i!Tnra  ( x  ,  y  ,  z  )  —  V  (  z  ) 


H3  ( x , y , z )  —4  I(z) 


(3.10) 


where  V(z)  and  I(z)  are  the  voltage  and  current  along  the  trans¬ 
mission  line,  i  is  the  propagation  constant  of  the  line,  k  is  its 
characteristic  impedance,  is  the  load  impedance,  Z^m  is  the 


characteristic  impedance  of  the  TE  mode: 

nm 

17  jw 

„nm  '  o 


(3.11) 


and  all  other  quantities  are  as  defined  previously.  Since  the 
transmission  line  analog  to  the  terminated  waveguide  is  a  line 
driven  at  z  =  -d  by  voltage  V(-d)  and  current  I ( -d )  and  termin¬ 
ated  at  z  =  -d  -  C  by  a  load  impedance  Z^,  the  solutions  to  the 

2 

transmission  line  equations  for  this  problem  : 


V ( z )  =  V(-d) 


I(z)  =  I(-d) 


Z^cosh  [  i  (C+d+z  )  ]  +  k  sinh  [  •»  ( C+d+z  )  ] 
Z^cosh[iC]  +  k  sinhOC] 


k  cosh  [*»(  C+d  +  z  )  ]  +  Z^sinh  [  t  (  C+d+z  )  ] 
k  cosh[nC]  +  Z^sinht+C] 


(3.12) 


give  the  following  solutions  to  the  corresponding  waveguide  prob¬ 
lem  using  (3.10) : 

E3Tim(x.y,z)  =  E3nm(x,y,-d)  F"m(z) 


H2Tnin(x  ,  y  ,  z  )  =  H3lim(  x  ,  y  ,  -d  )  fJJ®(z) 


(3.13) 


where 


ES71"1  (  x  ,  y ,  -d )  =  E3Tim(x,y) 


H3Tim(x,y,-d)  =  H3l7m(x,y) 


(3.14) 


F™(z) 


F^m(z) 


rj2  cosh  [rnm (  C+d+z  )  ]  +  Z^msinh  [rnm(  C+d+z  )  ] 
rj2  cosh[rnmCJ  +  Z™  sinh[rnmC] 

z"mcosh  [rnm(C+d+z)  ]  +  rj2  sinh  [rnm(C+d+z) ] 

Znmcosh[rnmC]  +  rj2  sinh[rnmC] 
o 


(3.15) 


o.  A.  Schelkunoff,  Electromagnetic  Waves,  Van  Nostrand,  Prince¬ 
ton,  NJ  (1943). 


and  the  amplitude  factors  $nm,  $nm,  {nm,  and  (ntn  in  H3nm(x,y)  and 

x  y  x  y 

E3um( x , y )  are  determined  by  the  magnitude  of  the  mode  fields  at  z 
=  -d.  The  same  procedure  using  the  well  known  expression  for  the 

input  impedance  of  a  loaded  transmission  line  of  length  C: 


Z.  =  k 
1 


Z^cosh(ifC)  +  k  sinh(iC) 


k  cosMiC)  +  Z^sinh(iC) 


(3.16) 


gives  the  impedance  of  the  TE  mode  at  z  =  -d  looking  into  M3: 
*  nm 


/7  3  nm  =  Znm 
'  o 


„  ,  . _nra_  >  _nm  .  ,  . .nm. . 

q2  cosh(r  C)  +  Zq  sinh(r  C) 


Znmcosh (rnmC)  +  q2  sinh(rnmC) 
o 


(3.17) 


where  all  quantities  are  as  defined  previously.  Both  (3.13)  and 
(3.17)  will  be  used  in  the  following  sections  to  obtain  the 
f ields  in  M3 . 

Mode  functions  have  many  other  properties  useful  in  the 
analysis  of  waveguide  problems.  Among  these  is  the  fact  that  they 
can  be  used  to  expand  arbitrary  electromagnetic  fields  as  series 
summed  over  an  infinite  set  of  discrete  modes  of  a  particular 
type.  That  is,  an  arbitrary  electromagnetic  field  with  no  elec¬ 
tric  field  component  in  the  z  direction  can  be  represented  by  a 
series  composed  of  TE  mode  functions.  Similarly,  an  arbitrary 
field  with  no  Hz  component  can  be  expressed  as  a  series  composed 
exclusively  of  TH  mode  functions.  Since  all  TE  and  TH  mode  func¬ 
tions  are  trigonometric,  both  types  of  expansions  will  be  two 
dimensional  Fourier  series  over  the  waveguide  cross  section.  For 


example,  if  E(x,y)  is  an  arbitrary  electric  field  with  no  z 


component  defined  over  the  waveguide  cross  section,  then  it  can 
be  represented  as  follows; 


00  00 


:(x,y)  =  )  )  C  e  (x,y) 

£  ^  nm  nm 


(3.18) 


n=0  m=0 


where  the  Fourier  coefficients  c  are 

nm 


nm 


A/2  B/2 

I  I 

-A/2  -B/2 


E  ( x  ,  y )  •  e  ( x , y )  dx  dy  . 

nm 


(3.19) 


The  latter  can  be  written  in  terms  of  the  x  and  y  components  of 
E(x,y)  as  follows: 

.y  (3.20) 


c  —  c  +  c  * 
nm  nm  nm 


where 


A/2  B/2 


x  ,, 
c  =  K 

nm  nm 


I  I 


E  (x,y)  e™(x,y)  dx  dy 

X  X 


-A/2  -B/2 


cy  =  K 
nm  nm 


A/2  B/2 

I  I 

-A/2  -B/2 


(3.21  ) 


Ey ( x , y )  e™m(x,y)  dx  dy 


and  e  (x,y)  and  e  (x,y)  are  given  by  (3.3).  Similarly,  the 
x  y 

transverse  components  of  an  arbitrary  magnetic  field  can  be 

expanded  in  a  Fourier  series  using  hTim(x,y). 

If  the  conductivity  of  the  waveguide  walls  is  finite,  then 
TE  and  TH  modes  are  coupled  at  the  walls  and  an  expansion  of  an 
arbitrary  field  in  terms  of  either  TE  or  TH  modes  alone  is  not 
possible.  Thus,  the  preceding  is  not  valid  in  the  general  case. 
However,  if  the  walls  are  composed  of  a  good  conductor  as  assumed 


$ 
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here,  it  has  been  shown  ’  that  TE  and  TH  mode  expansions  are 

valid  approximations  and,  with  one  exception,  the  other  mode 

characteristics  are  the  same  as  those  of  the  perfectly  conducting 

waveguide.  Thus,  the  characteristic  impedance  of  the  TE  mode  is 

nm 

still  given  by  (3.11)  and  the  cutoff  frequency  by  (3.7  ).  How¬ 
ever,  the  propagation  constant  rnm  is  modified  by  the  addition  of 
an  absorption  term  anm  that  depends  on  the  conductivity  of  the 

walls.  That  is,  (3.6)  is  replaced  by 

_  r  _  o  1/2  .  _ 


where 

nm 


nm . 2  .  .  .2 

>c  )  +  (J«) 


B  '/of1  +  (wcm/Jw)2 


(3.22) 


(3.23) 


.  mu  ,  . 

+  ( w  /iw 

c 


2 

2.21 

n  AB 

+ 

m  A 

n  B 

+ 

2.2 
m  A 

/ ,  .  B i  i  nm  .  .  . 2 

-  (1  +  j-)  (u>c  /jw) 


R  =  [w  p2/<2  o2)]l/i 

and  p2  and  o2  are  the  permeability  and  conductivity  of  the  wave¬ 


guide  walls.  When  o2  — ►  «>,  a' 


0,  and  (3.22)  reduces  to  (3.6). 


R.  E.  Collin,  Field  Theory  of  Guided  Waves,  McGraw-Hill,  New 
York  (1960). 
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V.  M.  Papadoppoulos ,  Quart.  J.  Mech.  and  Appl.Math.,  vol.7, 

(  1954  )  p  325 . 
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In  the  following  sections,  we  use  mode  functions  to  repre¬ 
sent  the  fields  first  at  the  M2/M3  interface  and  then  in  the  re¬ 
mainder  of  M3.  In  M3,  the  mode  fields  take  the  form  of  (3.7)  to 
account  for  reflections  between  the  termination  at  z  =  -(d  +  C) 
z  =  -(d  +  C)  and  the  driver  (incident  field)  at  z  =  -d. 

3 . 2  Mode  Expansions  for  Fields  Inside  a  Continuous 
Enclosure  Exposed  to  Arbitrary,  Time  Harmonic, 
External  Fields. 

The  first  step  in  computing  the  fields  in  M3  is  to 
expand  the  fields  incident  at  the  M2/M3  interface  in  terms  of 
mode  functions.  Since  the  fields  in  M2  are  TE  to  a  good  approxi¬ 
mation  and  since  the  dominante  mode  in  M3  is  also  TE ,  the  obvious 
choice  for  these  expansions  are  the  TE  mode  functions.  Expanding 

(2.29)  using  (3.2),  (3.3),  (3.4),  (3.18),  and  (3.21)  gives; 

00  00 

E2x(x,y,-d)  =  -  2  rj2  exp(-i2  d)  ^  ^  Cnmer(x,y* 

n=0  m= 1 
00  00 

H2y(x,y,-d)  =  2  exp(--r2  d)  ^  ^  cnmexB|x'y* 

n=0  m= 1 


(3.24) 


E2y ( x , y , -d )  =  2  q2  exp{-i2  d)  ^  ^  Cnmeym*x‘y* 

n=l  m=0 
00  00 

H2x(x,y,-d)  =  2  exp(-72  d)  ^  ^  cymeym(x,y) 


n=l  m=0 
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where 


A/2  B/2 

[  |  Hl®(x,y,0)  e^m(x,y)  dx  dy 


-A/2  -B/2 


(3.25) 


-  v2 


A/2  B/2 

I  I  H 

-A/2  -B/2 


c'  =  K  HI 

nm  nm  x 


Hl®(xfy,0)  e™(x,y)  dx  dy 

X  / 


and  all  other  quantities  are  as  defined  previously. 

In  this  representation,  each  field  in  (3.24)  is  the  sum  of 

an  infinite  series  of  plane  wave  fields  of  the  form 

E2nm(x,y,-d)  =  -2  q2  exp(--»2  d)  cX  enm(x,y) 
x  '  *  nm  x  . 


•  «  ry  HID  /  j 

H2y  ( x , y , -d 


2  exp(-i2  d)  cX  enm(x,y) 
nm  x 


(3.26) 


E2nm(x,y,-d)  =  2  n2  exp(-i2  d)  cy  enm  ( x  ,  y ) 

y  '  nmy 

H2nra( x , y , -d )  =  2  exp(-i2  d)  cy  enm(x,y) 

x  r  nmy 

for  0  <  n  <  »  and  0  <  m  <  «.  Each  of  these  fields  generates  re¬ 
flected  fields  in  M2  traveling  in  the  +  z  direction  and  transmit¬ 
ted  fields  in  M3  traveling  initially  in  the  -z  direction  and 
later  in  both  directions  after  reflections  between  the  back  and 
front  walls.  The  transmitted  mode  fields  in  M3  at  z  =  -d  gener¬ 
ated  by  each  incident  field  in  (3.26)  are  given  by 


(3.27) 


E3™(x,y,-d) 

=  E2nm 

X 

(x,y, 

-d) 

Tnm 

E 

H3^m (x,y,-d) 

=  H2nm 
y 

(x,y, 

-d) 

Tnm 

H 

E3"m(x,y,-d) 

=  E2nm 
y 

(  x  ,  y  , 

-d) 

Tnm 

E 

H3"m(x,y,-d) 

=  H2rUn 

X 

(  x  ,  y  , 

-d) 

rpHin 

h 

where  T^m  and  T„m  are  the  transmission  coefficients  for  the  elec- 
is  H 
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trie  and  magnetic  fields  of  the  TE  mode 

nm 

T"m  =  2  /7 3 nm/  (q2  +  q2™) 


2  q2/(q2  +  17 3  ) 


(3.28) 


q2  is  the  characteristic  impedance  of  the  wall  as  defined  previ¬ 
ously  and  n3nm  is  the  impedance  of  the  TE  mode  at  z  =  -d  look- 

nm 

ing  into  M3  as  given  by  (3.17). 

With  E3nm(x,y,-d)  and  H3iim  (  x  ,  y  ,  -d )  given  by  (3.27),  the 
mode  fields  at  all  other  locations  in  M3  can  be  obtained  from 
(3.13).  In  component  form,  the  result  is 
E3  ( x , y , z )  =  E3v  (x,y,-d)  F-  (z) 

X  XU 

iT^nni  /  v  IIOnm .  \ 

H3  ( x , y , z )  =  H3  (x,y,-d)  Fu  (z) 

y  y  n 

(3.29) 

_0nm.  .  _„nmy  „nm .  . 

E3y  ( x , y , z )  =  E3y  (x.y.-d)  F£  (z) 

H3x  (x,y,z)  =  H3x  ( x , y , -d )  FH  (z) 

where  the  fields  at  z  =-d  are 


E3nm 
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e"m(x,y) 

H3nm 

y 

(x,y,-d) 
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nm ,  . 

ey  ( x , y ) 

H3nm 
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(  x , y , -d ) 

nm 

5  X 

nm „  . 

ey  ( x , y ) 

the  amplitude  factors  are 


(3.30) 


.  hid  00  /  «  j  \  x  mHid 

$  =  -  2  n2  exp ( -•»  2  d)  c  T_ 

x  '  *  nm  E 


.  run  n  t  n  1  1  x  ^nm 
{  =  2  exp(-^2  d)  c  T.. 

y  nm  H 


?nm  =  2  q2  exp  ( -*»  2  d)  cy  t£® 
y  if  nm  E 


(3.31  ) 


.  nm  n  /  o  j  1  y  m n ® 

{  =2  exp  -t2  d)  c  T„ 

x  K  nm  H 

F^!m(z)  and  Ff)m(z)  are  given  by  (3.15),  and  all  other  quantities 
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are  as  defined  previously. 

Since  (3.29)  gives  the  fields  in  M3  due  to  the  nm-th  term 
in  the  series  expansion  of  (2.29),  the  total  field  at  any  point 
in  M3  due  to  fields  from  the  external  source  incident  at  the 
M2/M3  interface  is  equal  to  the  sum  of  terms  computed  from  (3.29) 
for  all  TE  modes  in  the  expansion  of  (2.29).  That  is, 


$x ( x , y , z )  =  ^  ^  E3^m(x,y,z) 


n=0  m=l 


( x  ,  y  ,  z  )  =  ^  S  H3"m(x,y,z) 


n=0  m=l 


(3.32) 


E3y(x,y,z)  =  ^  ^  E3^ra(x,y,z) 


n= 1  m=0 


H3x(x,t,z)  =  ^  ^  H3xm(x,y,z)  . 


n=l  m=0 

Since  (3.32)  satisfies  both  Maxwell’s  equations  at  all  points  in 
M3  and  all  boundary  conditions  at  the  M2/M3  interface,  it  repre¬ 
sents  a  unique  solution  to  the  shielding  problem  for  a  continuous 
rectangular  enclosure. 

Combining  (3.3),  (3.6),  (3.7),  (3.11),  (3.15),  (3.17), 

(3.22),  (3.23),  (3.25),  (3.28),  (3.29),  (3.30),  and  (3.31)  in 

(3.32)  gives  the  following  result: 

The  fields  at  any  point  inside  the  rectangular  enclosure 
shown  in  Figure  6  due  to  an  arbitrary  external  electromagnetic 
source  illuminating  a  single  wall  of  the  enclosure  can  be  written 


in  terms  of  the  tangential  components  of  the  source  magnetic 

field,  Hl®(x,y,0)  and  Hlfr(x,y,0),  incident  on  the  outside  surface 

— — — - ■  x -  y  — — — — - 

of  that  wall  as  follows: 


00  00 


E3x(x,y,z)  =  -  2  q2  exp(-t2  d)  ^  ^  Cnmexm(x,y)  TEm  FEm*z) 

n=0  m=l 
00  00 

H3y(x,y,z)  =  2  exp (-7  2  d)  ^  J  cV^^c.y)  T™  F^m(z) 


n  =  0  m=l 


OO  00 


(3.33) 


E3y(x,y,z)  =  2/72  exp(-i2  d)  S  \  c^n^y"1  ( x  '  y  >  T£m  FEm<z) 

n=l  m=0 

00  OO 

H3x(x,y,z)  =  2  exp ( -7  2  d)  ^  ^  Cnmeym(x,y)  FHm(z) 


n=l  m  =  0 


where 


q2  =  (ju  p2/o2) 


1/2 


72  =  (ju  p2  a2)1/2 
A/2  B/2 


x  „2 
c  =  K 
nm  nm 


cy  =  K2 
nm  nm 


Hl®(x,y,0)  e^®(x,y)  dx  dy 
y  x 


-A/2  -B/2 
A/2  B/2 


Hl®(x,y,0)  e™”(x,y)  dx  dy 
x  y 


K 


nm 


on 


e”m(x,y)  =  (mn/B)  cos [ nn ( x+A/2 ) /A ]  s i n [ ran ( y+B/2 ) /B ] 

nm  /  > 

e  (  x  ,  y )  =  - 


-A/2  -B/2 
(nn/A)2  +  (mn/B)^ 


=  1  when  n  =  0 
=  2  when  n  >  0 


-1/2 


[€  e  /(AB) ] 
on  on 


1/2 
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(nn/A)  s in [ nn ( x+A/2 ) /A )  cos [nm  ( y+B/2 ) /B ] 


ifj 


T"  =  2  r\2/ (rj2  +  q*  ) 

rj2  cosh(rnmC)  +  Znmsinh  (rnmC ) 
.nm  _nm  '  o 

no  -  z  ■  1  ■  ■  ~  .  i .  —  .i—i-i 

o  7nm  ,  .  nm. ,  „  ,  , _nm„ , 

Zq  cosh(r  C)  +  /7 2  sinh(r  C) 

q  ju 

.nm  'o 


nm  f,  nm.2  A  ..  .2 11'  . -1  J  nm 

r  =  (wc  )  +  ( jw)  (Co)  +  a 

w”m  =  7icf(n/A)2  +  (m/B)2! 


rj2  cosh  [rnm  (C+d+z )  ]  +  Znms inh [rnm ( C+d+z ) ] 

Fnm(z)  =  _ - _ 

E  _  ,  f_nm_.  .  7nm  .  ,  ,  .nm., 

q2  cosh[r  C]  +  Zq  sinh[r  C] 

Znmcosh [rnm( C+d+z ) ]  +  q 2  sinh [rnm(C+d+z ) ] 

Fnm(z)  =  _ ' _ 

H  „nm  ,  r,nm.,  ,  .  .  ,  r_nm_. 

ZQ  CGdh[r  C]  +  q 2  sinhtr  C] 

and  all  other  quantities  are  as  defined  previously. 


Mode  Expansions  for  Fields  Inside  a  Slotted 


Enclosure 


If  a  rectangular  slot  as  shown  in  figure  7  is  cut 
through  the  front  wall  of  the  enclosure  in  figure  6  and  the  rec¬ 
tangular  volume  of  the  slot  is  filled  with  a  material  of  perme¬ 
ability  /72’and  complex  permittivity  t2’,  then  the  latter  is 
transformed  from  a  continuous  to  a  discontinuous  enclosure,  If 
the  conductivity  o2’  of  the  new  material  is  significantly  less 


than  that  of  the  replaced  material  o2,  then  the  fields  reaching  a 
specified  point  in  M3  after  passing  through  the  slot  will  be  much 


larger  than  the  fields  at  the  same  point  before  the  slot  was  cut, 
provided  all  other  factors  are  unchanged.  If  the  coordinates  of 


the  center  of  the  slot  (x  ,y  )  and  its  dimensions  ( a,p )  are 

c  c 

known,  then  the  fields  in  M3  due  to  the  presence  of  the  slot  can 
be  computed  in  terms  of  the  products  of  the  tangential  compnents 
of  the  source  magnetic  field  incident  at  z  =  -d  and  the  form 
factors  of  the  slot.  The  method  used  is  exactly  the  same  as  in 
the  preceding  section  except  that  the  first  step  is  to  expand  the 
slot  fields  incident  at  z  =  -  d  (equation  (2.32)  rather  than 
equation  (2.29).  The  final  result  can  be  summarized  as  follows: 

The  fields  at  any  point  inside  the  rectangular  enclosure 
shown  in  figure  6  due  to  an  arbitrary  electromagnetic  source 
illuminating  a  rectangular  slot  located  in  the  front  wall  of  the 
enclosure  as  shown  in  figure  7  can  be  written  in  terms  of  the 
tangential  components  of  the  source  magnetic  field,  Hl^ ( x , y , 0 ) 

g 

and  HI  (x,y,0),  incident  on  the  slot  as  follows: 


E3x(x,y,z) 


00  00 

2  q2  exp(-i2  d)  S  S  ceX  enm(x,y) 
' x  ^  x  £  /  nm  x 

n  =  0  m=l 

00  00 


Tnm 

E 

x 


„nm 

FE 


(  z ) 


H3y(x,y,z) 


2  exp(-i2xd) 


V  V  ,  x  nm,  .  „nm  pnm,  , 

2  2  °hnm  (x'y)  TH  FH  Ul 

n  =  0  m=  1 

(3.35) 
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2  2  c  '!•«*> 
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H3x(x,y,z) 


=  2 


exp(-i2yd) 


ch’ 


nm 


nm ,  v 

ey  ( x , y ) 


nm 
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f-(z: 


n=l  m=0 
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where  q2  and  qQ  are  the  slot  impedances  when  it  is  driven  by  x 
directed  and  y  directed  electric  fields  respectively,  and  t2 

are  the  propagation  constants  of  the  fields  in  the  slot 
,2X  =  ju  u2/nz 


El 


T2y  =  jw  (j2/ q2  , 

and  fj2  are  the  effective  permeabilities  of  the  slot, 


ce 


nm 


ch 


nm 


ce 


nm 


chy 


nm 


nm 


nm 


K 


nm 


nm 


B/2  A/2 

Hly(x,y)  f ex (x , y , -d )  e^m(x,y)  dx  dy 
-B/2  -A/2 

B/2  A/2 

Hl®(x,y)  f h  ( x , y , -d )  e”m(x,y)  dx  dy 
y  y  X 

-B/2  -A/2 

B/2  A/2 

j  Hl®(x,y)  fey(x,y,-d)  eym(x,y)  dx  dy 
-B/2  -A/2 

B/2  A/2 

Hl®(x,y)  fhx(x,y,-d)  eym(x,y)  dx  dy 
-B/2  -A/2 


(3.36) 


(3.37) 


9 


f  e„  (x  ,  y ,  -d )  ,  f  hTr  (  x  ,  y ,  -d  )  ,  fe,r(x,y,-d) 

—A - J  - > - 

factors  as  defined  in  section  2.3, 


,nm 

E  = 

2 

q3nm/  ( rj2  +  q3nm) 

X 

,nm 

H  * 

X 

2 

172  /(n2  +  n3nm) 

,nm 

E  = 

2 

*7  3  /  (q2  +  773  ) 

y 

<7 

nm 

H  = 

2 

*72/(i7g  +  rj3nm)  , 

and  fh  (x,y)  are  form 
- x - - - 


(3.38) 


J*  V  V  -• 


_X^Ji  *  *  w  '  -  *  .- 


and  all  other  quantities  are  as  defined  in  the  preceding  section. 


If  there  are  two  or  more  slots  on  the  front  face  with  known 
dimensions  and  locations  and  if  the  impedances,  permeabilities, 
and  form  factors  of  these  slots  can  be  estimated,  then  the  fields 
at  a  given  location  in  M3  due  to  these  slots  can  be  obtained  by 
using  the  preceding  expressions  to  compute  the  fields  from  each 
slot  and  then  adding  the  contributions  from  all  slots  at  that 
location . 


4.  FIELDS  INSIDE  ENCLOSURES  EXPOSED  TO  SPATIALLY  UNIFORM 
EXTERNAL  FIELDS  GENERATED  BY  TIME  HARMONIC  SOURCES 

In  this  section,  the  general  expressions  for  the  fields  in¬ 
side  continuous  and  slotted  enclosures  are  evaluated  for  the  two 
simplest  cases  of  interest:  a  continuous  enclosure  exposed  to 
uniform  fields  over  its  front  surface,  and  a  slotted  enclosure 
exposed  to  uniform  fields  over  the  area  of  the  slot.  These  cases 
can  be  described  more  precisely  as  follows: 

Hl®<x,y,0)  a  HJ®(ju) 

s  s  (4*X) 

Hly ( x , y , 0 )  a 

over  the  following  areas: 

Continuous  Enclosure  (fig.  6) 

-  A/2  <  x  <  A/2 

-  B/2  <  y  <  B/2  (4.2) 

z  =  0 

Slotted  Enclosure  (fig.  6  and  7) 

x  -  a/2  <  x  <  x  +  a/2 
c  c 

yc  -  P/2  <  y  <  yQ  +  p/2  (4.3) 

z  =  0 

s  s 

where  HI  (jw)  and  HI  (jw)  are  functions  of  the  frequency  w  but 
x  y 

not  of  the  spatial  variables  x  and  y. 

These  cases  occur  typically  when  the  source  is  located  far 
from  the  enclosure  or  when  the  characteristic  dimension  of  the 
source  is  larger  than  that  of  the  enclosure  and  also  larger  than 
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the  separation  distance.  However,  they  can  also  occur  in  less  ob¬ 
vious  situations.  For  example,  it  was  shown  by  Monroe*  that  a 

g 

small  square  loop  antenna  10.6  in.  on  a  side  operating  at  10  Hz 
can  produce  a  nearly  uniform  magnetic  field  over  a  significant 
plane  area  oriented  perpendicular  to  the  plane  of  the  antenna  and 
located  1  m  from  its  center.  Furthermore,  since  the  dimensions  of 
most  slots  are  much  smaller  than  than  those  of  the  enclosure,  it 
is  clear  that  (4.1)  will  be  satisfied  in  many  cases  for  slotted 
enclosures  (4.3)  when  it  would  not  satsfied  for  continuous  enclo¬ 
sures  (4.2)  of  the  same  size  exposed  to  the  same  source.  Thus, 
(4.1),  (4.2),  and  (4.3)  can  describe  many  situations  of  practical 

importance . 


4.1  Continuous  Enclosure 


If  H1S  satisfies  (4.1),  then  E3  (x,y,z)  and 
x  y 

H3x(x,y,z)  must  be  independent  of  x  or  y  or  both  x  and  y.  But  it 

is  clear  from  (3.33)  and  (3.34)  that  in  order  for  E3  and  H3  to 

y  x 

be  independent  of  x  alone  or  both  x  and  y  it  is  necessary  that  n 

=  0,  which  implies  E3  =  H3  s  0.  The  latter  is,  of  course,  the 

y  x 

trivial  solution  that  will  only  occur  when  HI®  =  0.  Consequently, 
a  nontrivial  solution  requires  E3y  and  H3x  to  be  independent  of 

*R.  L.  Monroe,  A  Theory  of  Electromagnetic  Shielding  with 
Applications  to  MIL-STD  285,  IEEE-299,  and  EMP  Simulation, 

Harry  Diamond  Laboratories,  HDL-CR-85-052- 1 ,  Adelphi ,  MD 
(February  1985). 
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y.  This  implies  m  =  0.  Thus,  the  solution  in  this  case  is  obtain¬ 
ed  by  setting  in  :  0  in  the  expressions  for  E3  and  H3  in  (3.33). 

y  x 

g 

Similar  considerations  show  that  when  Hly  satisfies  (4.1) 

the  fields  E3  and  H3  are  independent  of  x  and  the  appropriate 
x  y 

expressions  for  these  quantities  are  obtained  by  setting  n  =  0 
in  (3.33). 

With  these  substitutions,  (3.33)  and  (3.34)  become 

00 

E3x(y,z)  =  -  2  q2  exp(--»2  d)  Hly(jw)  ^  c0mexm*y)  Te™  FEm(z) 


m=l 


H3y(y,z)  =  2  exp  ( -i  2  d)  Hly(jw)  ^  c0mexm(y>  FHm(z* 

m=  1 


(4.4) 


E3y(x,z)  =  2  rj2  exp(-»2  d)  Hl®(jw)  ^  cy0ey°(x)  tJ?0f£0<z) 


^nO^nO , 


n=  1 
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eyB(x)  =  sin [ nn ( x+A/2 ) /A ] 
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•v., 
-A 
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-1 


AVI 


and  all  other  quantities  are  either  as  previously  defined  or  are 
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While  the  expressions  in  (4.6)  remain  moderately  complicat¬ 
ed  in  spite  of  the  simplification  afforded  by  our  assumption  of 
spatially  uniform  incident  fields,  these  expressions  are  readily 
ammenable  to  evaluation  by  computer.  Calculations  based  on  (4.6) 
indicate  that  the  series  converge  rapidly  and  fewer  than  10  terms 
are  required  to  reduce  the  error  below  10%  at  most  locations  in 
the  enclosure.  In  general,  more  terms  are  required  at  locations 
close  to  the  side  walls  than  near  the  center  line  (x  =  y  =  0) 
where  five  terms  will  usually  suffice  for  10  %  accuracy.  Figures 
8  and  9  are  computer  plots  of  the  magnitude  of  H3^(y,z)  versus 
frequency  f  =  w/2n  at  locations  close  to  the  center  of  the  front 
wall  (x  =  y  =  0,  z  =  -0.05  m),  close  to  the  center  of  the  back 
wall  (x  =  y  =  0,  z  =  -6.05  m),  at  the  geometrical  center  of  the 
enclosure  (x  =  y  =  0,  z  =  -3.05  m),  and  near  the  center  of  the 

top  wall  (x  =  0,  y  =  1.75  m)  where  the  enclosure  parameters  are: 

-4 

A  =  6.1  m,  B  =  3.66  m,  C  =  6.1  m,  d  =  3.8  x  10  m, 

6  —  7 

o2  =  4.33  x  10  mhos/m,  fj 2  =  12.6  x  10  h/m, 


and  the  source  magnetic  field  is  harmonic  with  unit  magnitude: 


Hly(jw)  =  1.0  exp(jwt)  =  1.0  exp(j2nft) 

These  figures  show  that  the  general  trend  at 

ever  more  rapidly  decreasing  magnitude  with 

- 1/2 

reflecting  a  change  from  a  f  variation  to 

to  exp(-kf*^ 


(A/m)  . 

all  locations  is  an 
increasing  frequency 
f  *,  and  from  f-* 


).  Superimposed  on  this  trend  at  frequencies  above 
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Figure  8.  Magnitude  of  the  magnetic  field  inside 
a  continuous  enclosure  versus  frequency  at  locations 
near  the  center  of  the  front  (2=-. 05  m)  and  back 
walls  (z=-6.05  m). 


10  Hz  are  much  more  rapid  variations  corresponding  to  the  reson¬ 
ant  and  antiresonant  frequencies  of  the  enclosure.  Noteworthy 

7 

among  the  latter  is  the  pronounced  anti  resonance  at  f  a  2  x  10  Hz 
seen  near  the  center  of  the  front  wall  (x  =  y  =  0,  z  =  -0.05  m) 
but  nowhere  else.  Figure  8  indicates  a  spatial  variation  of  just 
under  one  order  of  magnitude  between  the  front  and  back  walls, 
and  figure  9  shows  a  variation  of  just  over  one  order  of  magni¬ 
tude  between  the  center  of  the  enclosure  and  a  point  (x  =  0,  y  = 
1.75  m,  z  =  -3.05  m)  0.08  m  from  the  center  of  the  side  wall.  The 
former  is  a  reflection  of  the  hyperbolic  variation  of  the  inter¬ 
nal  field  with  z,  and  the  latter  is  a  reflection  of  the  sinusoid¬ 
al  variation  of  the  internal  field  with  y  satisfying  the  boundary 
condition  H3  =  0  at  y  =  +1.83  m. 

y 

4.2  Enclosure  with  an  Inductive  Slot 

If  the  material  filling  the  slot  in  figure  7  is  free 
space  and  the  dimensions  of  the  slot  are  small  compared  to  the 
wavelength  A  of  the  fields  in  Ml,  then  the  permeability  of  the 
slot  is  equal  to  the  permeability  of  free  space  fj  and  the  slot 
impedances  are  pure  inductances.  That  is,  for  air-filled  slots 
where 

a  and  p  <  <  A  (4.7) 

the  slot  permeabilities  are  given  by 


and  the  slot  impedances  are 


1%  -  •>"Lx 


n 2  =  jwL 

y  y 


(4.9) 


where  Lx  is  the  inductance  of  the  slot  when  driven  by  an  electric 
field  in  the  x  direction  and  L  is  the  inductance  of  the  slot 

y 

when  driven  by  an  electric  field  in  the  y  direction.  The  propaga¬ 
tion  constants  for  the  fields  in  the  slot  are  derived  from  (4.8) 
and  (4.9)  as  follows: 


,2x  =  j»  fj2/rj2x  =  t-0/Lx 

,2y  =  J“  "F'72y  =  VLy 


(4.10) 


Thus,  the  general  expressions  for  the  fields  inside  a  slotted 

enclosure  can  be  specialized  to  the  case  of  an  inductive  slot  by 

substituting  (4.9)  and  (4.10)  into  (3.35). 

If  the  dimensions  a  and  /3  satisfy  the  relation  p  <<  a,  then 

figure  7  describes  a  narrow  rectangular  slot  with  its  long  axis 

parallel  to  the  x  axis  For  this  "horizontal”  air-filled  slot, 

the  inductances  satisfy  the  relation  L  <<  L  ,  and  similarly 

x  y 


rj2  <<  qZ 
x  y 


iZ  <<  72 
y  x 


(4.11) 


by  virtue  of  (4.9)  and  (4.10).  Applying  (4.11)  to  (3.35)  gives 


|E3x  J  <<  |E3y| 

lH3yl  <<  l«3x| 


(4.12) 


provided  |H1^|  a  |Hly|.  That  is,  E3y  and  H3x  are  the  dominant 
field  components  inside  an  enclosure  with  a  narrow  horizontal 


slot  provided  j H 1 S  |  is  not  very  much  larger  than  |H1S|. 

y  x 

Similarly,  if  the  slot  dimensions  satisfy  a  <<  p,  then  fig¬ 
ure  7  describes  a  narrow  vertical  slot  with  its  long  axis  paral¬ 
lel  to  the  y  axis.  In  this  case,  the  slot  inductances  satisfy 

L  <<  L  and  (4.11)  becomes 
y  x 

<<  nl 

(4.13) 

i2  <<  t2 

x  y 

The  dominant  fields  inside  the  enclosure  are  now  E3  and  H3  . 

x  y 

That  is, 

lE3yl  <<  lE3*l 

(4.14) 

|H3XI  <<  |H3yl 

provided  |H1*|  *  j  HI ®  |  . 

The  preceding  shows  that,  in  general,  the  fields  inside  an 
enclosure  with  a  single  narrow,  air-filled,  horizontal  or  verti¬ 
cal  slot  will  be  dominated  by  a  single  electric  field  component 
and  a  single  magnetic  field  component  and  that  the  remaining  com¬ 
ponents  will  be  negligible  by  comparison.  In  the  following,  we 
conrider  the  horizontal  slot  and  obtain  explicit  expressions  for 

E3  and  H3  .  The  reader  will  then  be  able  to  transform  these  into 
y  x 

the  corresponding  expressions  for  E3^  and  H3y  inside  an  enclosure 
with  a  vertical  slot. 

Substituting  for  q2  from  (4.9)  and  for  from  (4.10)  into 

(3.35)  and  separating  the  m  =  0  terms  in  the  sum,  we  obtain 


E3y(x,y,z)  =  2juLy  exp ( -pQd/Ly )  ^  cen  ey(x)  TE  FE ( Z ' 


,  V  V  y  nm .  v  _nm  _nm ,  k 

+  2  2  nm  ey  <x'yl  TE  FE  <z) 


n=  1  m= 1 


(4.15) 


(x,y,z)  =  2  exp(-pod/Ly)  ^  chy  ey(x)  Tjj  fJJ(z) 

n=  1  y 

00  00 

,  Y  Y  i_y  nm.  v  -nm  rnm.  v 
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where 


n=l  m= 1 


B/2  A/2 


:ey  =  K2  HlS(jw) 
nm  nm  x 


chy  =  K2  HlS(jw) 
nm  nm  x 


-B/2 

-A/2 

B/2 

A/2 

-B/2 

-A/2 

fey(x,y,-d)  ey  (x,y)  dx  dy 


fhy(x,y,-d)  ey  (x,y)  dx  dy 


e  (  x  ,  y )  =  -  (  n7i/ / A  )  s  in  [  nrr  (  x+A/2  ) /A  ]  cos  [nm  (  y+B/2  )  /B  ] 

y 

T™  =  2  /73nm/(jwL  +  rj3nm)  (4.16) 

y  y 

tH®  =  2  jwL  /  ( jwL  +  rj3nm) 
n  y  y  y 

y  y 

ce  =  ce  „ 

n  nm  m=0 

chy  =  chy 

n  nm  m  =  0 

n ,  .  nm .  , 

ey(x)  =  ey  ( x , y )  m  =  0 

nO  nm  m  =  0 

and  all  other  quantities  are  defined  in  section  3. 

To  determine  E3  and  H3  from  the  preceding  expressions, 

y  x 

it  is  necessary  first  to  specify  the  form  factors  fe  (x,y,-d)  and 


V 


fhy(x,y,-d)  which  describe  the  spatial  variation  of  the  incident 


fields  at  z 


-d  for  -  A/2  <  x  <  A/2  and  -B/2  <  y  <  B/2  and 


then  to  use  these  factors  to  evaluate  the  Fourier  coefficients 

y  y 

ce  and  ch  by  carrying  out  the  indicated  integrations.  How- 
nm  nm  , 

ever,  exact  expressions  for  fe  and  fh  are  unknown,  and  there- 

y  y 

fore  in  order  to  proceed,  we  must  resort  to  approximations.  In 
making  these  approximations,  we  will  be  guided  by  the  fact  that 
the  incident  fields  at  z  =  -d  are  determined  by  the  fields  inside 
the  slot  (xc  -  a/2  <  x  <  xc  +  a/2,  yc  -  p/2  <  y  <  yQ  +  p/2) 
and  by  the  boundary  conditions  over  the  remainder  of  the  z  =  -d 
surface.  Since  the  slot,  like  the  enclosure,  is  equivalent  to  a 
section  of  rectangular  wavegude  where  TE  fields  dominate,  it 
follows  from  the  results  quoted  in  section  3  that  the  spatial 
variation  of  the  slot  fields  can  be  written  as  sums  using  TE  mode 
functions  appropriate  to  the  slot.  That  is, 

00  oo 


e2yU,y)  ^  "sr  ''y^*'5'’ 


8=1  r=0 


00  00 


(4.17) 


(X,y)  =  ^  Jfsr  *yr<x’y> 


s=l  r=0 

where  xy  and  *y  are  the  Fourier  coefficients  of  the  expansion 
sr  sr 


k  (x,y)  =  v  sinfsnfx-x  +  a/2 )/a]  cos[rn(y-y  +  p/2)/p)  (4.18) 

y  c  c 


with 


u  -  -  7i  s  /  a  . 

However,  there  are  no  propagating  modes  in  this  expansion  due  to 
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our  assumption  that  the  dimensions  of  the  slot  are  small  compared 
to  the  wavelength  of  the  source  field,  and  consequently  the  field 
in  the  slot  will  be  determined  by  relatively  few  low-order  TE 
modes  with  the  TE^q  mode  playing  a  dominant  role.  Thus,  a  first 
order  approximation  to  the  structure  of  the  fields  in  the  slot 


can  be  obtained  from  (4.17): 


10 


e2  (x,y)  =  h2  ( x , y )  s  k  (x,y) 
y  ^  y 


=  sin[n(x-x  +a/2)/a]  . 
c 


(4.19) 


Of  course,  this  is  just  equivalent  to  approximating  e2  and  h2 

y  x 


by  the  first  terms  in  their  mode  expansions.  It  implies  that  the 
fields  in  the  slot  have  the  structure  of  a  half-sine  wave  in  the 
x  direction  and  are  uniform  in  the  y  direction.  This  appears  to 
be  a  reasonable  assumption  for  the  fields  in  a  narrow  slot  excit¬ 
ed  by  a  uniform  external  field.  Thus  the  structure  of  the  inci¬ 
dent  field  at  z  =  -d  can  be  approximated  by  combining  the  struc¬ 
ture  of  the  slot  field  from  (4.19)  with  continuous  extensions  of 
the  incident  slot  field  from  the  slot  to  the  remainder  of  the  z  = 
-d  plane.  In  the  case  of  the  electric  field,  there  are  no  contin¬ 
uous  extensions  beyond  the  vicinity  of  the  slot  since  the  remain¬ 
der  of  the  z  =  -d  plane  is  composed  of  a  good  conductor.  Thus, 


fey  can  be  approximated  by 


fe  (x,y,-d)  =  p(x)  q  (y) 
y  C 


(4.20) 


where 


sin[n(x  -  x  +  a/2)/a]  for  x  -  a/2  <  x  <  x  +  a/2 
c  c  "  ■  c 


for  -A/2  <  x  <  x  -  a/2  and  x  +  a/2  <  x  <  A/2 


P(x) 


1 


q^(y)  = 


for  y  -  p/2  <  y  <  y  +  p/2 
c  c 

0  for  -B/2  <  y  <  yc  -  p/2  and  yc  +  p/2  <  y  <  B/2 


In  the  case  of  the  magnetic  field,  there  are  continuous  exten¬ 
sions  of  the  incident  field  beyond  the  slot  in  the  +y  directions 
since  the  slot  field  at  (x,+/5/2,-d)  will  drive  currents  in  these 
directions  on  the  z  =  -d  plane.  These  currents  will  be  reflected 
at  the  side  walls  with  a  reflection  coefficient  of  nearly  +1  re¬ 
sulting  in  an  increase  in  H2x(x,y,-d)  in  the  vicinity  of  y  =  + 

B/2.  A  precise  determination  of  fhy(x,y,-d)  would  model  the  z  =  - 
d  plane  in  the  y  direction  as  an  asymmetrically  driven  transmis¬ 
sion  line  with  appropriate  loads  at  y  =  +  B/2.  Here,  however,  we 


will  approximate  fhy  analytically  as  follows: 
f hy ( x , y , -d )  =  p ( x )  qh(y) 


(4.21) 


where  p(x)  is  the  same  sinusoidal  function  used  in  (4.19),  q^(y) 


is  defined  by 

■  s,y  +  q,  for  -  B/2  <  y  <  yc  -  p/2 
qh(y)  =  1  for  yc  -  p/2  <  y  <  yc  +  p/2 

8,y  +  q2  for  +  p/2  <  y  <  B  /2 


s,  =  -  l/(yc  -  p/2  +  B/2) 
s,  =  -  l/(yc  +  p/2  -  B/2) 


q i  =  1  +  ( yc  -  P! 2)/(yc  -  p/2  +  B/2) 
q,  =  1  +  (yc  +  p/ 2)/(yc  +  p/2  -  B/2) 


(4.22) 


(4.23) 
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Figure  10  shows  plots  of  p(x)  and  q^(y)  for  A  =  10,  B  =  8,  a  =  6, 

p  =0  .5,  xc  =  +1  and  yc  =  -2  where  the  units  are  arbitrary. 

When  (4.20)  and  (4.21)  are  substituted  into  (4.16),  the 

integrals  can  be  evaluated  in  closed  form  to  give  explicit  ex- 

y  y 

pressions  for  the  Fourier  coefficients  ce  and  ch  .  These  can 

nm  nm 

then  be  combined  with  (4.15)  in  the  following  statement: 

The  principal  components  of  the  field  at  any  point  inside 


the  enclosure  shown  in  figure  6  due  to  a  uniform  electromagnetic 


source  illuminating  a  narrow,  horizontal,  air-filled  slot  with 


center  coordinates  x  and  y  and  dimensions  a  and  p,  as  shown  in 


figure  7,  can  written  in  terms  of  the  horizontal  component  of  the 


source  magnetic  field  Hlx(jw)  incident  on  the  slot  as  follows: 


E3y(x,y,z)  =  2 jwLy  exp( -/uQd/Ly )  ^  cej  ey(x)  te  Fe(z) 

n=l  y 

l  1  c 


n=l  m=l 


(4.24) 


( x , y , z )  =  2  exp ( -pQd/L  )  \  chy  e"(x)  T^  F^(z) 


00  00 


V  V  ,  y  nm-  -  -.nn  „nni-  > 

2  2  nm  (x’yl  TH  FH  ,z> 


n= 1  m=l  J 

where  L  is  the  slot  inductance  when  driven  by  a  y-directed 

- y - I - - - 

electric  field 

cey  =  -  HlS(jw)  p  P  K2  nn/A 
n  x  ^  n  n 

chy  =  -  HlS(jw)  Q  P  K2  nw/A  I 

n  x  o  n  n 


(4.25) 


■  -  *  >  •  -  '  V'  ■  •  V  «  ‘  •  •  •  V.-W  V  V.V  V  S'-V.VA 


s  *. 


I  I  I 

-5  -4 


i  i  i 

-3  -2 


Figure  10.  Plots  of  p(x)  and  q^(y)  for  A  =  10. 

B  :  8,  a  :  6,  p  :  .5,  x  =  1 ,  and  y  =  -2. 

c  c 

(Arbitrary  units). 
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HI  ju  P  Qe  K  nn/A 

x  n  m  nm 

HlS(jw)  P  Qh  K2  nn/A 
x  n  m  nm 


sinfn(A-na) (x  +a/2)/(aA)] 


a/2 )  /a 


sin[n(A+na  (x  +a/2)/(aA)j  -  n(x  -a/2 )/a  +  nn/2] 
c  c 

A  +  na 


(4.26) 


sin [ n ( A-na )( x  -a/ 2)/(aA)]  -  n(x  -a/2)/a  -  nn/2] 
c  c 

A-na 


sin[n( A+na) (x  -a/2)/(aA)]  -  n(x  -a/2)/a  +  nn/2 
c  c 


A  +  na 


provided  A  *  na 


Qo  =  p  +  (yc  -  p/2  +  B/2 )  s1(yc-  p/2  -  B/2)/2  +  qj 
-  (yc  +  p/2  -  B/2)  s2(yc  +  p/2  +  B/2)/2  + 


(4.27) 


sin[mn(y  +  p/2  +  B/2)/B]  -  sin[mn(y  -  p/2  +  B/2)/B] 


(4.28) 


Qh  =  Qe  + 
m  m 


[s,(yc  -  p/2)  +  q  4 ]  sin[mn(yc  -  p/2  +  B/2)/B] 
mn/B 


Sjfcoslnmfy  -  p/2  +  B/2)/B]  -  1] 

+  - E _ 

(nu/B) 2 

(4.29) 

[s, (yc  +  p/2)  +  q  2 ]  sin[mn(yc  +  p/2  +  B/2)/B] 
mn/B 

Sj[cos(m7i)  -  cos[mn(y  +  p/2  +  B/2)/B] 

+  - E _ 

( nm/B )  2 

8i»  s2»  Qi i  and  q2  are  given  by  (4.23)  and  all  other  quantities 
are  as  defined  previously. 

Although  the  preceding  are  much  more  complicated  than  the 

corresponding  expressions  for  a  continuous  enclosure,  they  are 

just  as  amenable  to  evaluation  as  the  latter  with  the  aid  of  a 

computer.  Figures  11  and  12  are  computer  plots  of  |H3x(x,y,z)| 

versus  frequency  at  several  locations  inside  an  enclosure  with  a 

narrow  horizontal  slot  at  the  center  of  one  wall  (x  =  y  =  0) 

c  'c 

exposed  to  a  uniform  harmonic  source  of  unit  magnitude 

(exp(  j27ift )  )  .  The  enclosure  and  slot  parameters  used  in  these 

plots  are 

A  =  B  =  C  =  2.7  m,  d  =  6.4xlO-3  m 

o2  =  4.33X10^  mhos/m,  p2  =  12.6x10  ^  h/m 

a  =0.5  m,  p  =  1.6X10-3  m,  Ly  =  3.2xl0-8  h 
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Figure  11.  Magnitude  of  the  magnetic  field  inside 
an  enclosure  with  a  narrow  air-filled  slot  at 
points  near  the  center  of  the  slot  (z  =  -  .05  m) 
and  the  center  of  the  back  wall  (z  =  -  2.65  m). 
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Figure  12.  Magnitude  of  the  magnetic  field  inside  an 
enclosure  wi*h  a  narrow  air-filled  slot  at  a  point 
near  the  center  of  one  side  (x  =  1.30  m)  and  at  the 
center  of  the  enclosure  (x  =  0.00  m). 
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which  is  derived  from  figure  7  of  Monroe  (1973)  .  A  total  of  156 

terms  were  used  in  the  calculations  by  truncating  the  series  at  n 
=  m  =  12.  Plots  obtained  using  more  terms  were  found  to  be  indis¬ 
tinguishable  from  the  ones  in  figures  11  and  12.  Plots  using 
fewer  than  156  terms  were  found  to  differ  from  those  in  figures 
11  and  12  by  scale  factors  that  are  only  weakly  dependent  on 
frequency.  For  example,  the  lower  of  the  two  closely  spaced 
curves  in  figure  13  is  a  replot  of  the  curve  labeled  z  =  -  0.05  m 
in  figure  11  where  156  terms  were  used  consisting  of  144  terms  in 
the  double  sum  over  n  and  m  plus  12  terms  from  the  single  sum 
over  n.  This  curve  is  nearly  identical  in  form  to  the  upper  curve 
which  was  obtained  using  only  the  12  terms  in  the  sum  over  n. 

The  figures  suggest  that  magnetic  fields  at  all  locations 
are  virtually  independent  of  frequency  below  the  lowest  cutoff 
frequency  of  the  slotted  enclosure.  This  contrasts  with  the  con¬ 
tinuous  enclosure  where  the  field  is  a  monotonical ly  decreasing 
function  of  frequency  below  cutoff.  The  fields  in  the  slotted  en¬ 
closure  also  show  a  significantly  greater  variation  with  position 
than  was  seen  in  the  continuous  enclosure.  Thus, in  figure  11,  the 
front  to  back  change  is  two  orders  of  magnitude  compared  to  just 
under  one  order  of  magnitude  change  in  figure  8  in  spite  of  the 
fact  that  the  slotted  enclosure  is  much  smaller.  Evidently  fields 

*^R.  L.  Monroe,  EMP  Shielding  Effectiveness  and  MIL-STD  285, 
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REQUEN 
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i 


Pi 


in  the  slotted  enclosure,  besides  being  much  larger  in  magnitude, 
are  more  localized  than  in  the  continuous  enclosure. 

Above  the  cutoff  frequency,  the  magnetic  field  shows  basic¬ 
ally  the  same  resonant  behavior  as  in  the  continuous  enclosure 
but  with  vastly  larger  magnitudes.  In  fact,  the  figures  indicate 
that  the  interior  field  can  exceed  the  exterior  field  at  resonant 
frequencies.  This  suggests  that  a  0.5  m  Blot  makes  this  enclosure 
worse  than  useless  as  an  electromagnetic  shield  at  these  frequen- 


5.  TIME  DOMAIN  MAGNETIC  FIELDS  INSIDE  ENCLOSURES  DUE  TO 
SPATIALLY  UNIFORM,  EXTERNAL,  TRANSIENT  SOURCES 


If  an  external  source  generates 
components  that  are  known  functions  hi 
t  with  LaPlace  transforms: 

OD 


Hl®(j»)  =  L[hl®(t ) ] 


exp(jut)  hi 


transient  magnetic  field 

s(t)  and  hls(t)  of  the  time 
x  y 

s(t)  dt 
x 


0 


Hl®(jw)  =  L[hl*(t)] 


r00 

exp(jut)  hl^r(t)  dt 


(5.1) 


J0 

that  are  spatially  uniform  over  the  areas  defined  by  (4.2)  for 
continuous  enclosures  and  by  (4.3)  for  slotted  enclosures,  then 
expressions  for  the  transient  fields  inside  continuous  and  slot¬ 
ted  enclosures  exposed  to  this  source  can  be  obtained  by  substi¬ 
tuting  (5.1)  into  (4.6)  and  (4.24)  and  computing  the  inverse  La- 
Place  transform  of  the  resulting  expression.  Since  (4.6)  and 
(4.24)  are  infinite  series  each  term  of  which  is  a  function  of  ju 
,  this  process  will  consist  of  a  terra-by-term  application  of  the 
inverse  LaPlace  transform  to  the  series.  If  the  inverse  transform 
of  each  term  in  the  series  can  be  evaluated,  then  the  result  will 
be  new  series  each  term  of  which  is  a  function  of  t  as  well  as 
the  internal  spatial  variables. 


In  the  following  sections,  we  use  this  procedure  to  compute 
magnetic  fields  inside  enclosures  exposed  to  transient  fields  in 


the  form  of  decaying  exponentials,  damped  sinusoids,  and  ration¬ 
alized  exponentials-  a  group  of  functions  that  alone  or  in  com¬ 
bination  can  be  used  to  approximate  the  waveforms  of  most  natural 
and  man-made  sources.  Although  a  rigorous  mathematical  justifica¬ 
tion  of  term-by-term  inversion  of  infinite  series  is  «  difficult 

1 5 

and  important  problem  as  Doetsch  has  emphasized,  we  will  not 
attempt  such  a  justification  here.  Our  justification  will  rest  on 
the  fact  that  the  time  domain  series  show  rapid  convergence  and 
produce  waveforms  that  are  physically  plausible  and  consistent 
with  measurements. 

5.1  Continuous  Enclosure 


If  H3  (x,z)  and  H3  (y,z)  from  (4.6)  are  rewritten  in 
x  y 


the  following  form 

OO 

H3x(x,y;jw)  =  ^Sx(x,z;ju) 


n=  1 


(5.2) 


H3y(y,z;jw)  =  2Sy(y,z:jw) 


m=  1 


where  the  frequency  is  displayed  explicitly  as  an  independent 


variable  through  the  factor  jw  : 


s”  (  x  ,  z  ;  jw  )  =  2exp[-r1/2(  jw)  1/2]  H 1  ®  (  jw  )  cVJ  (  x  )  tJJ  (  jw  )  Fjj(z;jw) 


(5.3) 
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Sy  <  y ,  z  ;  jw )  =  2exp[-r1/2( jw)1/2]  HI®  (  jw  )c*e“(  y  )tJJ(  jw  )  F®(z:jw) 


r  =  do2/;2  =  4T 


(T  =  characteristic  diffusion  time)  (5.4) 


then  the  inverse  LaPlace  transforms  of  H3x(x,z;jw)  and 
H3y(y,z;jw)  are 


h3x(x,z;t)  =  L  1 [H3x(x,z; jw) ]  =  ^h3x(x,z;t) 

n=  1 
00 

h3  (y,z;t)  =  L-1[H3  (y,z;jw)]  =  Sh3®(y,z;t) 

y  y  y 


(5.5) 


where 

h3x(x,z;t)  =  L_1[Sx(x,z; ju) ] 

h3y(y,z;t)  =  L~ 1 [ Sy ( y , z ; jw ) ] 
and  the  inverse  transform  is  defined  as  follows 


(5.6) 


L  l  G  (  x  f  y ,  z  ;  s  )  ]  = 


expfts)  G(x,y,z;s)  ds 


(5.7) 


where  jw  is  replaced  by  the  complex  variable  s  =  u  +  jv. 

Thus  the  problem  of  computing  time  domain  magnetic  fields  inside 

this  enclosure  reduces  to  an  evaluation  of  a  sequence  of  inverse 

LaPlace  transforms  (5.6)  where  Sn(x,z;jw)  and  Sm(y,z;jw)  are 

x  y 

given  by  (5.3). 

Since  S^(x,z;jw)  and  Sn(y,z;jw)  are  both  complicated  functions 
x  y 

of  jw ,  for  all  n  and  m,  we  shall  not  attempt  an  exact  evaluation 
of  these  transforms.  Instead,  we  shall  use  approximations  to  re¬ 
duce  Sn(x,z; jw)  and  Sm(y,z;jw)  to  the  product  of  a  pair  of 
x  y 
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functions  each  of  which  has  an  exact  transform.  This  will  allow 


us  to  express  h3n(x,z;t)  and  h3m(y,z;t)  as  sequences  of  convolu- 

x  y 

tion  integrals  in  the  time  domain  which  can  readily  be  computed 

to  give  h3  (x,z;t)  and  h3  (y,z;t). 
x  y 

The  approximation  is  based  on  the  fact  that  the  exponential 
factor  in  (5.3)  effectively  prevents  frequencies  above  a  critical 
value  Fc  =  1  /  (  tt  r  )  from  reaching  the  interior  of  the  enclosure*. 
Since  calculations  show  that  F  is  well  below  the  lowest  cutoff 


c 


frequencies  for 


Fc  =  1  / ( nr  ) 


it  follows  that 

Ji±L  <<  i 


n 


all  room  sized  shielded  enclosures 

(  f  “c  '<2">ln=1  =  V(2A) 

(  l  “c  /<2">Ib=1  =  C0/(2B) 
only  frequencies  satisfying 


(5.8) 


and 


(5.9) 


u 


m 


<<  1 


can  affect  the 

_m  n  ,  _m  „ 
Z  ,  r  ,  and  r  f 
o 

proximations : 


„n  .  n 

Z  -  n  jw/w 
o  '  o  c 

rj  ni  .  m 

Z  =17  Ju/« 
o  o  c 


interior 
rom  (4.6) 


magnetic  fields  in  this  case.  Thus,  Z^, 
can  be  replaced  by  the  following  ap- 


*R.  L.  Monroe,  A  Theory  of  Electromagnetic  Shielding  with 
Applications  to  MIL-STD  285,  IEEE-299,  and  EMP  Simulation, 
Harry  Diamond  Laboratories,  HDL-CR-8 5 -052- 1 ,  Adelphi ,  MD 


(  February  1985)  . 


(5.10) 


rn  =  »n/c 

c  o 


m  m 
r  =  w  /C 
c  o 


where 
n 


w 


c 

m 


jjC  n/A 
o 


(5.11) 


u  =  nC  m/B 
c  o 


and  the  absorption  terms  in  rn  and  rm  have  been  neglected.  With 


these  approximations,  calculations  show  that 


r)3n  2  z”  tanh(rnC) 


/7 3m  2  Z®  tanh(rmC) 


(5.12) 


and 


Tjj(jw)  » 


2(p2/o2) 


1/2 


(fj2/o2)1/Z  +  ( jw)  1/2(/7q/w")  tanh[rnC] 


„m , 


2 ( p2/o2 ) 


1/2 


T“(j«)  a 


(p2/o2)1//2  +  ( jw)  1/2(i7o/w“)  tanh[rmC] 


F^(z;jw)  =  cosh  [ /"n  (  C+d+z  )] /cosh  [rnC] 


(5.13) 


F®(z;jw)  =  cosh[/"m  (C+d+z)  ] /cosh[rnC]  . 


Substituting  (5.13)  into  (5.3)  and  the  latter  into  (5.6),  we  ob¬ 
tain 


h3"(x,z;t)  =  L  1 [ R j ( x , z ; jw )  Rj(jw)3 


h3™ ( y , z ; t )  =  L_1(R™(y,z; jw)  R®(jw>] 


(5.14) 


where 


R*?(x,z;jw)  =  2  HlS(jw)  cyen(x)  cosh  [  fn  (  C  +  d  +  z  )) /cosh  [rnC  ] 


2(^2/a2)l/£‘  exp[-r<  jw)*' ‘  ] 


1/  £. 


RgUw)  = 


(jj2/a2  )  +  (  jw  )  ly^2  (/Jq/w”  )  tanh[rnC] 


(5.15) 


R. (y.z; jw)  =  2  HI  (jw)  c  e  (y)  cosh[r  (C+d+z ) ] /cosh[r  C] 
l  y  m  x 


R^Jw)  = 


2(^2/o2)1/2  exp[-r( jw) 1/2] 


(p2/o2 ) 1^2  +  ( jw ) 1 ^2 ( 77  /w® )  tanh [rmC] 


Since  L  * [ R^ ( x , z ; jw ) ]  and  L  *  [ R™  (  y  ,  z  ;  jw  )  ]  can  be  written  in  terms 


of  the  incident  time  domain  source  fields  (5.1) 


L_1[H1®( jw) ] 


L_1[L[hl®(t) ] ] 


h  1 S  ( t ) 
x 


and 


L-1[Hl®(jw)]  =  L-1[L[hl®(t) ] ]  =  hl®(t) 

J  J  J 


which  are  assumed  to  be  known  functions  of  t  for  t  >  0  and  since 


L  ^[R^jw)]  and  L  ^[RgUw)]  are  also  known  [Campbell  and  Foster, 


No.  809] 


16 


L  1  [  R  g  (  j  )  ]  =  2  rn(t) 


L  ^ [ R™ ( jw ) ]  =  2  rm(t) 


(5.16) 


where 


n 


t’1/2  r2  1 

r  (t)  =  1/2  n  exp('  n1  *  ■  h-2 

7Tb  (  b  ) 


exp 


n 


(bn ) 2 


1 G 

G.  A.  Campbell  and  R.M.  Foster,  Fourier  Integrals  for 
Practical  Applications,  D.  Van  Nostrand  Company,  Inc. 


Princeton,  N.J.  (1948). 


.  ■«.  \  %  %  ». 


%  -%  •-  "u  *.  V  V  *1. 


x  erfc 


r  ^  t 

IT17*  1 


(5.17) 


rm(t) 


-1/2 


exp(-  ^-)  - 


r  t1/2 

X  iT177  +  "b5" 


n  tanh[rC] 
n  'o 


w“ ( p2/o2 ) 


f7Qtanh  [  r  C  ] 


(5.18) 


w"'(p2/o2) 


and  erfc[f(t)]  is  the  complementary  error  function 


erfc[f(t)]  s  —\~j2  exp(-v2)  dv 
ti  1 ' 

f  (t) 


(5.19) 


it  follows  from  the  convolution  theorem  [Doetsch,  No.  23] 15  that 

the  inverse  transforms  in  (5.4)  can  be  written  in  terms  of 

hl®(t),  hl®(t),  rn ( t ) ,  and  rm(t).  That  is, 
x  y 


h3"(x,z;t)  =  4  cyen(x)  In(t)  cosh [ rn ( C+d+z )] /cosh ( rnC ] 
x  n  y 

h3™ ( y , z ; t )  =  4  c*e“(y)  Ira(t)  cosh [ rm ( C+d+z ) J /cosh [ rmC ] 


(5.20) 


G.  Doetsch,  Guide  to  the  Application  of  LaPlace  Transforms,  D, 
Van  Nostrand  and  Company  Ltd.,  London  (1961). 


92 


where  the  convolution  integrals  I  (t)  and  I  (t)  are 
ft 

In(t)  =  hl°(t-g)  rn(g)  dg 


Im(t)  = 


fhly(t 

Jo 


(5.21) 


-g)  r  (g)  dg  . 


Substituting  (5.20)  into  (5.5)  gives  the  desired  expressions  for 

h3  (x,z;t)  and  h3  (y,z;t) .  The  complete  statement  is  as  follows: 
x  y 

Approximate  expressions  for  the  time  domain  magnetic  fields 
at  any  point  inside  a  continuous  enclosure  due  to  a  transient 
external  electromagnetic  source  generating  spatially  uniform 
fields  on  one  wall  of  the  enclosure  can  be  written  in  terras  of 

g 

the  tangential  components  of  the  source  magnetic  field  hl^ ( t )  and 

g 

hi  (t)  incident  on  the  outside  surface  of  that  wall  as  follows 

- y - 

oo 

h3  (x,z;t)  =  4  \  cyen(x)  T  n  ( t )  cosh  [  r11  (  C+d+z  )] /cosh  [  rnC  ) 
x  Z.  n  y 

n=  1 


(5.22) 


(y,z;t)  =  4  S  cXem(y)  I m ( t )  cosh [ rm ( C+d +z )] /cosh [ rmC ] 
L  m  x 


where 


In(t)  = 

•  0 

Im(t)  = 

•  0 


hi  (t-g)  r  (g)  dg 


hl„(t-g)  r  ( g )  dg 


r  ( g )  and  r  (g)  are  given  by  (5.17),  T  and  r  by  (5.10),  and  the 
remaining  quantities  are  as  defined  in  section  4.1. 


>  .*  s'  s^S's  r,V.‘  s'  V_  w  s  sju 


jv  t>t>’ *>'  *  rv.v.v..  v:  *7 


h 


i 

- 

; 

» 

i 

* 


I 

i 


i 

i 

i 


i 


i 


s  s 

When  hi  (t)  and  hi  (t)  are  specified,  the  convolution  int- 
x  y 

egrals  can  be  evaluated  for  each  n  and  m,  and  h3  (x,z;t)  and 
h3y(y,z;t)  can  be  computed  using  a  truncated  version  of  (5.22). 
However,  the  integrals  cannot  usually  be  reduced  to  closed-form 
expressions,  and,  therefore,  must  be  evaluated  using  a  numerical 
integration  technique.  Since  these  techniques  replace  integrals 
with  sums  over  the  time  interval  of  interest  (0,t)  for  each  term 
in  the  series,  the  single  sums  in  (5.22)  become  double  sums. 
Calculations  based  on  (5.22)  where  the  integrals  have  been  eval¬ 
uated  using  the  trapezoidal  rule  with  40  points  per  interval  and 
the  series  are  truncated  at  n  =  m  =  10  show  excellent  convergence 
for  most  locations  and  time  intervals.  This  approach  was  used  to 
obtain  the  waveforms  in  figures  14,  15,  16,  and  17,  where  the  in¬ 

ternal  fields  h 3  (solid  curves)  are  multiplied  by  a  scale  factor 

y 

4  s 

of  10  and  the  incident  fields  hl^  (dotted  curves)  are  given  by 
the  following: 

Double  exponential  hlS(t)  =  exp(-4xl0^t)  -  exp(-4xl0®t) 

sy  9  7 

Rationalized  exponential  hi  (t)  =  l/[exp(-3xl0  )  +  exp(2.3xl0  t)] 

Damped  sinusoid  (1.59  MHz)  hlS(t)  =  exp(-4xl0^t.)  sin(lxlO^t) 

C  n 

Damped  sinusoid  (6.37  MHz)  hl^(t)  =  exp(-2.5xl0  )  sin(4xl()  t) 
the  enclosure  is  characterized  by 

A  =  6.1  m  B  =  3.66  m  C  =  6.1  m  d  =  3.8X10  ^  m 

o2  =  4.33x10^  mhos/m  p/2  =  12.6x10  ^  h/m 
and  the  location  is  the  center  of  the  enclosure 
x  =  y  =  0.00  m  ,  z  -  -  3.05  m 
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Figure  15.  Transient  magnetic  field  h3y  x  10  A/»  (solid 

curve)  at  the  center  of  a  continuous  enclosure  with  one 


wall  exposed  to  hl“  (dotted  curve)  with  a  rationalized 
exponential  waveform. 
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INTERNAL  MAGNETIC  FIELB  (h3v)  IN  AHPS/M  TIMES  10(4) 


FOR  ENCLOSURE  PARAMETERS  AND  INCIDENT  FIELD  SHOWN  BEL 
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Figure.  17  Transient  magnetic  field  h3  x  10  A/m  (Bolid 

y 

curve)  at  the  center  of  a  continuous  enclosure  with  one  wall 
exposed  to  hls  (dotted  curve)  with  an  exponentially  damped, 

y 

1.59  MHz  sinusoidal  waveform. 


These  waveforms  agree  with  earlier  results  of  Monroe1  who 
based  his  calculations  on  a  one-term  approximation  to  the  trans¬ 


ient  field  at  z  =  -d.  They  show  that  T^,  the  time  required  for 
h3y  to  reach  its  maximum,  is  greater  than  the  characteristic  dif¬ 


fusion  time  of  the  enclosure  T  ,  where 

c 

Tc  s  r/4  =  a 2  p2  d2/4 


(5.23) 


2  200  ns, 


(1  ns  =1x10  s) 


for  all  four  incident  waveforms  and  is  independent  of  T^,  the 

s 

time  reqiured  for  hl^.  to  reach  its  peak.  For  these  cases,  T^ 

g 

depends  primarily  on  the  duration  of  hl^  with  a  longer  pulse  pro¬ 
ducing  larger  value  of  T^,  that  is,  a  slower  "rise"  time.  This  is 
seen  most  clearly  in  figures  14  and  15  where  the  double-exponen¬ 
tial  pulse  with  a  duration  of  1000  ns  produces  an  internal  field 
with  Tg  s  820  ns  while  the  rationalized  exponential  with  a  dura¬ 
tion  of  200  ns  produces  T^  =  420  ns.  These  figures  also  confirm 
Monroe’s  calculations  showing  that  the  maximum  value  of  h3  de- 

y 

pends  strongly  on  the  pulse  duration  with  long  incident  pulses 
generating  much  larger  internal  fields  than  comparable  short  pul¬ 
ses.  Thus,  the  double  exponential  with  a  maximum  value  (1  A/m) 
equal  to  that  of  rationalized  exponential  generates  an  internal 
field  with  a  maximum  value  nearly  five  times  larger.  This  is  due 
to  the  fact  that  longer  pulses  carry  more  energy  at  low  frequen¬ 
cies  which,  as  shown  by  figures  8  and  9,  are  generally  less  at¬ 
tenuated  in  passing  through  the  wall  of  a  continuous  enclosure. 


R.  L.  Monroe,  A  Theory  of  Electromagnetic  Shielding  with 
Applications  to  MIL-STD  285,  IEEE-299,  and  EMP  Simulation, 
Harry  Diamond  Laboratories,  HDL-CR-85-052- 1 ,  Adelphi,  MD 
(February  1985). 


The  importance  of  the  frequency  content  of  the  incident  pulse  in 
determining  the  magnitude  of  the  internal  field  is  strikingly 
illustrated  by  figures  16  and  17  where  in  both  figures  the  inci¬ 
dent  pulse  is  an  exponentially  damped  sine  wave.  Here,  the  1.59 
MHz  sine  wave  (Figure  17)  generates  an  internal  field  three  times 
larger  than  the  field  generated  by  the  6.37  MHz  sine  wave  despite 
the  fact  that  its  amplitude  is  smaller  and  its  damping  greater. 

The  accuracy  of  the  approximations  used  to  obtain  (5.22) 
can  be  gauged  by  computing  | H 3  (y,z;ju)|  from  (5.2)  with  T™(jw) 

y  n 

and  F®(z:jw)  determined  by  (5.13)  and  comparing  the  result  with 
|H3  (y,z;jw)|  computed  with  the  original  expressions  (4.6).  This 
is  done  in  figure  18,  where  the  upper  curve  is  obtained  from 
(4.6)  and  the  lower  curve  is  the  approximation.  For  the  scale 
used  here,  the  two  curves  are  indistinguishable  at  frequencies 
below  20  MHz.  Above  20  MHz,  the  curves  diverge  with  the  approxi¬ 
mation  falling  farther  below  the  exact  curve  until  the  first  en¬ 
closure  resonance  at  approximately  70  MHz  is  reached.  At  frequen¬ 
cies  above  the  resonance,  the  curves  partially  converge  but  then 
diverge  again  in  the  vicinity  of  the  second  resonance.  This  al¬ 
ternating  divergence  and  convergence  occurs  at  every  resonance 
and  is  superimposed  on  a  general  downward  trend  that  is  virtually 
the  same  for  both  curves.  Clearly,  the  principal  difference  be¬ 
tween  the  exact  and  approximate  expressions  is  that  the  former 
accounts  for  the  effects  of  enclosure  resonances  while  the  latter 
does  not.  However,  for  the  case  shown  in  figure  18  and  for  most 
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ft  CONTINUOUS  ENCLOSURE  FOR  THE  ENCLOSURE 
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Figure  18.  Plots  of  | H3y ( x , y , z  ; j» ) |  computed  from  equation 
(4.6)  using  exact  expressions  for  T^(jw)  and  FH(z;jw)  (upper 
curve)  and  using  approximate  expressions  (5.13)  for  T^(jw) 
and  F®(z;jw)  (lower  curve). 


cases  of  interest,  the  lowest  enclosure  resonance  occurs  at  a 


frequency  where  the  field  is  highly  attenuated  with  respect  to 
the  external  field,  and  higher  frequency  fields  are  attenuated 
even  more.  Thus,  enclosure  resonances  are  likely  to  have  little 
effect  on  the  transient  internal  field  inside  most  continuous 
enclosures,  and  the  fact  that  the  approximation  fails  to  account 
for  resonances  will  have  little  effect  on  the  accuracy  of  the 
transient  responses  computed  with  this  approximation.  Significant 
errors  may  occur  in  the  case  of  a  relatively  large,  thin-walled 
enclosure  where  condition  (5.8)  is  violated.  When  (5.8)  is  vio¬ 
lated,  the  field  at  resonance  will  be  less  attenuated  than  the 
fields  at  most  lower  frequencies,  and  the  transient  reponse  of 
the  enclosure  may  be  dominated  by  energy  coupled  into  the  higher 
frequencies.  This  is  most  likely  to  occur  when  the  incident  field 
is  a  narrow-banded  pulse  centered  on,  or  near,  a  resonant  fre¬ 
quency.  In  such  a  case,  one  should  not  expect  the  approximation 
to  yield  an  accurate  transient  reponse. 

For  the  enclosure  used  to  obtain  figures  14,  15,  16,  and 

17,  the  critical  frequency  Fc  equals  0.4  MHz  and  Co/(2B)  =  41 
MHz.  Hence,  (5.8)  is  easily  satisfied,  and  the  approximation  used 


5.1  Slotted  Enclosure 


As  in  the  preceding  section,  we  begin  calculating  the 
transient  response  of  a  slotted  enclosure  by  rewriting  the  gener¬ 
al  frequency  domain  expression  for  the  internal  magnetic  field 
(4.24)  in  a  form  that  is  convenient  for  term-by-term  application 
of  the  inverse  Laplace  transform.  That  is,  we  rewrite  H3x(x,y,z) 
from  (4.24)  as  follows 

00  OO  00 

H3x(x,y,z; jw)  =  ^Sx(x,z;jw)  +  ^  ^  S™ ( x , y , z ; jw )  (5.24) 


n= 1  m=l 


where 


s"(x,z;jw)  =  2  Hl®(jw)  cy  e"(x)  tJJ  (jw)  f"(z;jw) 
x  x  n  y  n  n 

y 

n  /  <  i  o  tt  i  ^  \  y  nm  /  ,  Mnm  /  .  ,  ^ n m .  ,  ■ 

SY  (x,y,z;jw)  =  2  HI  (jw)  c "  e  (x,y)  T  (jw)  F  (z;jw) 
x  x  run  y  n  n 

c5  =  -  exp ( -p  d/L  )  Q  P  nn/A 
n  r  ro  y  o  n  n 

cy  =  -  exp(-p  d/L  )  P  Qh  nw/A 

nm  *  ro  y  n  m  nm 


(5.25) 


(5.26) 


and  apply  the  inverse  transform  to  this  expression  to  obtain  the 
transient  magnetic  field: 


h3x(x,y,z;t)  =  L  [ H3x ( x , y , z ; jw ) 1 


=  ^  h3x ( x , z ; t )  +  ^  ^  h3xm(x,y,z;t) 


(5.27) 


n=l  m=l 


where 


h3” ( x , z ; t )  =  L-1 [s"(x,z; jw) ] 
h3xm(x,y,z;t)  =  L" 1 [ S™* x , y , z ; jw ) ] 


(5.28) 


Thus,  to  compute  the  transient  internal  field,  it  is  necessary  to 
eva luctc  the  sequence  of  inverse  LaPlace  transforms  (5.28)  where 


S^(x,z;jw)  and  S^m ( x , y , z ; jw )  are  given  by  (5.25)  and  (5.26). 

As  in  the  case  of  the  continuous  enclosure,  it  will  be  nec¬ 
essary  to  introduce  approximations  in  order  to  reduce  S^(x,z;jw) 
and  S™m ( x , y , z ; jw )  to  expressions  that  will  allow  us  to  evaluate 
their  inverse  transforms  in  terms  of  convolution  integrals  for 
all  n  and  m.  Since  figures  11  and  12  show  that  enclosure  reson¬ 
ances  can  play  an  important  role  in  determining  the  largest  field 
inside  a  slotted  enclosure  over  a  wide  range  of  locations  and 
frequencies,  it  is  clear  that  approximations  (5.9)  cannot  be  ex¬ 
tended  to  the  slotted  enclosure.  Thus,  Zn ,  Z^m,  rn,  and  rnm  can- 

o  0 

not  be  replaced  by  approximations  based  on  (5.9).  Instead,  we 
must  use 

Zo  =  /7njw/?n(^w) 


„ni  .  .  nm,  .  . 

Zo  =  nQ jw/s  (jw) 


_  n  n  <  •  \  i  r>  i  ^ 

T  =  ?  (jw)/Cq  +  a 

nm  nm .  .  ,  nm 

r  =  $  ( jw  )  /C  +  a 

o 


(5.29 


where 


. n ,  .  .  r  ,  n .  2  .  .  .2.1/2 

$  ( JW )  =  [ ( wc )  +  ( JW )  ] 

nm  .  .  .  . j  nm .2  .  .  .2,1/2 

$  (jw)  =  [ (w  )  +  (jw)  ] 


w  =  7T C  n/A 
c  o 


wnm  =  nC  [(n/A)2  +  (m/B)2 
c  o 


(5.30) 


(5.31  ) 


v.v.v 


WPS* ss*** 


'.'.'S' 


>»a 


These  are  the  original  expressions  rewritten  in  a  form  that  emph¬ 
asizes  their  frequency  dependence  through  the  factors  $n(jw)  and 
?nm(jw).  With  Z",  Z™,  rn,  and  rnm  given  by  the  preceding,  the 
remaining  approximations  used  for  the  continuous  enclosure  can  be 
applied  to  the  slotted  enclosure  without  any  additional  loss  of 
accuracy.  Thus,  we  can  again  replace  r;3n,  /7 3 nm ,  F[|(z;jw),  and 
F^m(z;jw)  by 

/7 3 n  a  z"  tanh  [  rnC  ] 


„nm  „nm  ,  nm  , 

n 3  a  Z  tanh[r  C] 

o 

(5.32) 

F™(z;jw)  a  cosh [rn < C+d+z )] /cosh [fnC ] 

H 

F[!m(z;jw)  at  cosh  [  rnm(  C+d+z  )] /cosh  [  rnmC  ] 

n 

With  these  approximations  for  77 3 n  and  ^7 3nm ,  we  obtain 

2  jwL 

Tn(ju)  =  _ - _ 

H  j«Ly  +  z"tanh(rnC] 

2 jwL  (5'33) 

Ojw)  =  - Z _ 

H  jwL  +  Znmtanh[/-nmC] 

y  o 

from  (4.16).  Combining  (5.33)  with  F[l(z;jw)  and  F!]m(z;ju)  from 


(5.33) 


(5.32)  gives 


2jwL  coshtf  (C+d+z)] 

tJJu»>  f[|  (  z  ;  jw  )  =  - l — - - - - — 

jwLycosh [ rC ]  +  Z^sinh[rC] 

(5.34) 

2jwL  cosh [ rnm ( C+d+z ) ] 

-.nro .  .  ,  rnni ,  y 

T„  (jw)  F„  z;  jw  =  - - - 

H  H  v  .  .  ,  r  _  nm  _ .  .  „nm  .  ,  r.niD„. 

jwLycosh[r  C]  +  Zq  sinh[r  C] 

Substituting  (5.34)  into  (5.25)  and  the  latter  into  (5.28),  we 
obtain 

h3^(x,z;t)  =  L_1tRj(x;jw)  R^fzjjw)] 


**’  i'* ,*■  s  *  •  *  »V'*  a  *  a  *  " 


h3^n’(x,y,z;t)  =  L-1 [R^m(x,y; jw)  R^z'.jw)] 


(5.35) 


where 


R"(x;jw)  =  2  c^ey(x) 


R2(z;jw)  = 


2  L  Jn(ju)  cosh[$n(jw)  (C+d  +  z)/C  ] 

_ y _ o _ 

Ly?n(  jw  )cosh  [  ?n  (  jw  )C/Cq]  +  r7osinh[$n(  jw)C/CQ] 


n  m  .  ■  i  a  y  n  id  .  .  . .  -  s  ,  ,  . 

Rj  (X,y;jw)  =  2  c^me  (x,y)  Hlx<jw) 


(5.36) 


2  L  <nm(jw)  cosh[cnm( jw) (C+d+z)/C  ] 
y  o 


Dnm y  ' "  o 

R?  (  z  ;  jw  )  =  - - - - - 

"  Ly?nm(  jw)coshUnm(  jw)C/CQ]  +  /7osinh[?nm(  jw)C/Co) 

and  ,  Z^m,  rn,  and  rnra  have  been  replaced  by  (5.29)  and  (5.30) 
after  dropping  the  absorption  terms  a11  and  anm  from  the  latter. 

Following  the  same  approach  used  for  the  continuous  enclo¬ 
sure,  we  will  compute  the  inverse  transforms  of  R”(x;jw), 
R^  (  z  ;  jw  )  ,  R^m  (  x  ,  y  ;  jw  )  ,  and  R^fzjjw)  and  use  these  quantities  to 
write  h3x(x,z;t)  and  h 3^™ ( x , y , z ; jw )  in  terms  of  sequences  of  con¬ 
volution  integrals  by  applying  the  convolution  theorem  to  (5.35). 

s  s 

Since  hi  (t)  is  the  inverse  transform  of  HI  (jw),  we  obtain  the 
x  x 

inverse  transforms  of  R^(x;jw)  and  R^fx.yjjw)  immediately; 


r  j  (  x  ;  t )  =  L-1[Rj(x;  jw)  ]  =  2  c^ey(x)  hl^U) 


r?ra  ( x  ,  y  ;  t )  =  b~ 1  [  (  x  ,  y ;  jw  )  ]  =  2  cy  enm(x,y)  hls(t) 

i  l  nm  y  x 


(5.37) 


However,  the  inverse  transforms  of  R^(z;jw)  and  R^fzjjw)  are  not 
as  readily  available.  These  transforms  do  not  appear  in  any  of 
the  standard  collections,  and,  therefore,  they  must  be  calculated 


directly  without  recourse  to  tabulated  functions.  The  most 


.  V_v_  -  .V.  S_  ^ 


straightforward  way  to  do  this  is  by  carrying  out  partial  frac¬ 
tion  expansions  of  R^izjjw)  and  R^lziju)  and  applying  the  in¬ 
verse  transform  to  each  term  in  the  expansions.  Since  R^tzis)  and 
R^izjs)  are  meromorphic  functions^  of  the  complex  variable  s, 
these  expansions  and  their  inverse  transforms  will  be  infinite 
series.  The  derivation  in  appendix  A.  gives 


r>n,  .  ,  2  +  ,  n.2,,..2 

R2<z;jw)  =  -  -  ^  +  2  («c)  +  (jw) 


Wk(z, 


(v>c)  +  (u>k)  +  (ju) 


R£  ( z ; jw 


2  f  .  0[.  nm . 2  , 2 

y  ■+  i  +  2  («c  )  +  (jw) 


(5.38) 


Wk(z) 


k=i  <«c")  +  (V  +  <•>">' 


where  unm  and  w11  are  given  by  (5.31), 
c  c 


L  C 
y  o 


Vzl  = 


2  ip  cos[p.  (C+d+z)/C  ] 
k  o 

*  [cos(p,  )  -  p.sinfp.  )]  +  coslp,  ) 


(5.39) 


w.  =  p.  C  /C 
k  k  o 

and  p^  k  —  1,2,3*. 


are  roots  of  the  equation 


*  p  cos (p)  +  sin(p)  =  0 


(5.40) 


Applying  the  inverse  transform  to  (5.38)  and  using  transform 


E.  C.  Titchmarsh,  The  Theory  of  Functions,  2d  Ed.,  Oxford 
University  Press  (1939). 
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pairs  33  and  38  in  Doetsch  gives 

-1  r^n 


(  z  ;  t )  =  L_  *  [  Rg  (  z  ;  ju> )  ]  =  if  Oft)  -  2  ^ 


00  *“k*2  Wk^Z^ 


k=  1 


(5.41) 


r2m(  z  ;  t )  =  L-1 [R2m( z ; jw ) ]  =  t  o(t)  -  2  ^ 


®  /  1 2  i  >  j  nm ,  l 

“k  Wk  2  8in(w}t  t> 


k  =  1 


w, 


nm 


where  6(t)  is  the  delta  function  and 


,  n . 2  ,  .  .21 

(wc)  +  <«k> 


1/2 


w. 


nm 


,  nm  2  ,  .2 

(w  )  +  (wk) 


,  1/2 


(5.42) 


We  can  now  apply  the  convolution  theorem  to  (5.35)  using 
(5.37)  and  (5.40)  to  obtain 


h3n(x,z;t)  =  2  cyen(x)  In(z;t) 
x  n  y 


rnm , 


i  0nm .  .  r\  y  nm ,  »  _ nm .  ,  . 

h3x  (  x  ,  y ,  z  ;  t )  =  2  c^e  (x,y)  I  (z;t) 


(5.43) 


where 


In(z;t)  = 


Tnm .  ... 

I  (  z  ;  t )  = 


hl^(t-g)  r^fzjg)  dg 


,t 


(5.44) 


,  i s .  .  nm ,  _  x  , 

hlx(t-g)  r2  (z;g)  dg 


0 


These  expressions  can  now  be  substituted  into  (5.27)  to  determine 
h3  (x,y,z;t).  The  final  result  can  be  stated  as  follows: 


Doetsch,  Guide  to  the  Application  of  LaPlace  Transforms,  D. 


Van  Nostrand  and  Company  Ltd.,  London  (1961). 
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An  approximate  expression  for  the  time  domain  magnetic 


field  at  any  point  inside  an  enclosure  with  a  narrow,  horizontal, 


air-filled  slot,  as  shown  in  figure  7,  due  to  a  transient  exter¬ 


nal  source  generating  a  spatially  uniform  field  over  the  surface 


of  the  slot  can  be  written  in  terms  of  the  tangential  component 


of  the  source  magnetic  field  hi  (t)  incident  on  the  slot  as  fol 


h3x(x,y,z;t)  =  2  ^cye^(x)  I'*(z;t) 


(5.45 ) 


_  V  V  y  nm,  Tm,  .  . 

+  2  2  2  Cnmey  (x>y;  1  (z;t) 


where 


In(z;t) 


n=l  m=l 


h 1 x ( t-g )  r2(z;g)  dg 


Tnm  . 

I  ( z 


-t 

5t)  = 

•  o 


hl®(t-g)  r2m(z;g)  dg 


cy  and  cy  are  given  by  (5.26),  rl!|z;t!  and  r!!m  ar~  given  by 
— n - nm -  2 - 2 - 

(5.41),  and  e[^(x)  and  e„  ( x ,  y )  are  as  defined  previously. 


Since  r2(z;t)  and  r^  (z;t)  are  infinite  series  for  each  n 
and  m,  In(z;t)  and  Inm(z;t)  will  also  be  infinite  series  for  each 
n  and  m.  Consequently,  the  single  sum  in  (5.45)  is  actually  a 
double  sum,  and  the  double  sum  is  a  triple  sum.  When  the  func¬ 
tional  form  of  hlx<t)  i9  specified,  In(z*,t)  and  Inm(z;t)  can  be 
computed  and  h3^ ( x , y , z ; t )  can  be  written  explicitly  in  terras  of 
these  sums.  For  example,  if  hl®(t)  can  be  written  in  the  form  of 


the  difference  between  two  decaying  exponentials 

g 

hi  =  h  [exp(-^.t)  -  exp(-7„t)] 


(5.46) 


.  ( x ,  y ,  z ;  t )  =  2  ^  ^  c^e^(x)  l”(z;t) 


n=l  k=0 


00  CO  00 


(5.47) 


*  2  1  1  1  Cnmeym(X’y)  IC“<zit) 


n= 1  m=l  k=0 


where 


1 0  (  Z  1 1  >  =  Jo™*2’1*  =  *  hlx^^  =  *  hQ  [  exp  (  -7  j  t )  -  exp(-72t)] 
l"(z;t)  =  -  2(wk)2Wk(z)  j£<t)/wk 
I^m(z;t)  =  -  2(«k)2Wk(z)  Jkm(t)/wjm 

(5.48) 

n  wk  exP(-1,1t)  +  '»19in(«j|t)  -  w”  cos(u^t) 

Jk(t)  =  *  ho  -  - 

(V  +  (V 

<*>k  exp( -7  gt )  +  72sin(«kt)  -  w£  cos(u^t) 


1 2'  +  K> 


nm  ,  .  .  .  ,  nm .  v  nm  .  nnr  . 

w,  exp(-7,t)  +  71sin(w,_  t)  -  uk  co s(u,  t) 


/  v "  i  /  run . 

(,i'  *  <uk  » 

m  exp ( - 7 gt )  +  72sin(wk 


Similar  expressions  can  be  obtained  when  hl^(t)  assumes  other 
functional  forms,  such  as  the  damped  sinusoid.  However,  when  the 
convolution  integrals  cannot  be  evaluated  in  closed  form,  one 
must  resort  to  numerical  integration.  This  has  the  effect  of 
replacing  double  sums  with  triple  sums  and  triple  sums  with  quad¬ 


ruple  sums  and  can  greatly  complicate  the  process  of  computing 


h3x(x,y,z;t) . 


The  solid  curves  in  figures  19,  20,  21,  22,  23,  24,  25,  and 

26  are  plots  of  h3  (x,y,z,;t)  computed  with  (5.47)  and  (5.48)  for 

the  same  enclosure  used  to  obtain  figures  11,  12,  and  13  where 

the  center  of  the  slot  coincides  with  the  center  line  of  the 

enclosure  (x  =  y  =0).  These  curves  represent  time  histories  of 
c  c 

the  magnetic  field  at  points  on  the  center  line  of  the  enclosure 

ranging  from  the  center  of  the  back  surface  of  the  slot  (x  =  y  = 
_  3 

0,  z  =  -6.44x10  m)  to  the  center  of  the  back  wall  (x  =  y  =  0,  z 

=  -2.7  m).  The  incident  field  (5.46)  is  characterized  by  h  =  1.1 

o 

7-1  8—1 

A/m,  k  j  =  1  x  10  s  ,  and  t ^  -  5  x  10  s  .  These  parameters 

produce  a  pulse  (dotted  curve)  with  a  peak  field  of  1  A/m  at  t  = 

8  ns.  It  should  be  noted  that  h3x (x , y , z ; t )  and  hl®(t)  are  plotted 

to  the  same  vertical  scale  in  figures  19  and  20,  but  h3  (x,y,z;t) 

is  multiplied  by  a  factor  of  10  in  the  remaining  figures.  The 

_  9 

time  scale  is  in  nanoseconds  (ns)  where  1  ns  =  1  x  10  s,  and  t  = 
0  ns  corresponds  to  the  arrival  of  the  incident  field  at  z  =  0. 

The  computer  time  required  to  complete  the  calculations  for 
these  figures  was  significantly  reduced  by  including  only  the  sum 
over  n  and  k  in  (5.47),  that  is,  by  dropping  the  triple  sum  over 
n,  m,  and  k.  The  approximation  is  justified  in  this  case  by 
figure  13,  which  shows  that  the  difference  between  the  exact  and 
approximate  expressions  is  relatively  small  and  virtually  inde¬ 
pendent  of  frequency  at  a  location  (x  =  y  =  0,  z  =  -  0.05  m) 
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Figure  19.  Transient  magnetic  field  h3x  (A/m)  inside 

_  O 

a  slotted  enclosure  at  the  point  (x=y=0 , z=-6 . 4X10  m) 
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Figure  21.  Transient  magnetic  field  h3x  X  10  (A/m) 
inside  a  slotted  enclosure  at  the  point 
( x=y=0 , z=- . 338  m ) . 
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Figure  22.  Transient  magnetic  field  h3x  X  10  (A/m) 
inside  a  slotted  enclosure  at  the  point 
( x=y=0 ,  z  =  - . 675  m ) 
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Figure  23.  Transient  magnetic  field  h3x  X  10  (A/m) 
inside  a  slotted  enclosure  at  the  point 
(x=y=0,z=-1.35  m). 
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Figure  24.  Transient  magnetic  field  h3x  X  10  (A/m) 
inside  a  slotted  enclosure  at  the  point 
(x=y=0 , z=-2 . 03  m ) . 
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Figure  25.  Transient  magnetic  field  h3x  X  10  (A/m) 
inside  a  slotted  enclosure  at  the  point 
( x=y=0 , z=-2 . 7  m ) . 
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close  to  the  slot.  Moreover,  other  calculations  for  this  case 
show  that  the  difference  becomes  even  smaller  as  the  field  point 
moves  farther  away  from  the  slot.  Thus,  we  can  expect  figures  19 
to  26  to  accurately  represent  the  wave  form  of  the  transient 
field  while  somewhat  overestimating  its  magnitude  at  locations 
close  to  the  slot.  These  curves  were  generated  by  truncating  the 
series  at  n  =  12  and  k  =72  and  calculating  h3^  at  time  intervals 
of  .5  ns  using  a  PC-XT  operating  at  4.77  MHz.  A  70  ns  record  re¬ 
quired  a  running  time  of  approximately  11  hours. 

The  field  at  the  slot  as  shown  in  figure  19  is  consistent 
with  the  one  term  approximation  of  Monroe  *  in  that  it  consists 
of  a  unipolar  ( dc )  pulse  carrying  a  smaller  bipolar  { ac )  pulse 
that  recurs  with  certain  variations  at  a  fixed  interval  (18+  ns) 
corresponding  to  the  time  required  by  a  signal  moving  at  approxi¬ 
mately  the  speed  of  light  to  reflect  from  the  back  wall  and  re¬ 
turn  to  the  front.  Since  the  same  features  are  seen  at  other 
points  in  the  enclosure,  we  will  follow  Monroe  and  refer  to  these 
as  the  static  or  stationary  field  and  the  propagating  field.  As 
shown  in  figure  19,  the  form  of  the  static  field  is  close  to  that 
of  the  incident  field  and,  like  the  latter,  it  approaches  zero 
exponentially  as  t  -»  ».  On  the  other  hand,  the  propagating  field 
has  a  much  more  complicated  structure  that  changes  with  each 
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reflection  and  bears  no  obvious  resemblance  to  the  incident 
field.  Clearly,  the  static  field  is  derived  mainly  from  those 
frequencies  in  the  incident  pulse  that  are  below  the  lowest  TE 
cutoff  frequency  for  the  enclosure  (flat  portion  of  the  curves  in 
fig.  11),  while  the  propagating  field  is  derived  from  frequencies 
above  the  cutoff.  Thus  the  propagating  field  will  be  continually 
reflected  between  the  front  and  back  walls  and  will  continue  to 
propagate  after  the  static  field  has  decayed  to  insignificance. 

Since  the  static  field  is  composed  of  frequencies  below  TE 
cutoff,  we  would  expect  it  to  decay  rapidly  with  distance  from 
the  slot.  This  expectation  is  borne  out  by  figures  21,  22,  and 
23  where  it  is  seen  that  the  peak  static  field  at  z  =  -0.68  m  is 
only  l/20th  of  the  peak  field  at  the  surface  of  the  slot.  When 
the  field  point  moves  an  additional  0.65  m  from  the  slot  to  the 
center  of  the  enclosure  as  in  figure  24,  the  static  field  ib  re¬ 
duced  to  a  barely  visible  level.  Thus  the  static  field  is  large 
compared  to  the  propagating  field  only  at  points  very  close  to 
the  slot.  At  most  locations  in  the  enclosure,  the  propagating 
field  dominates  as  shown  by  figures  24,  25  and  26. 

The  rapid  decay  of  the  static  field  with  distance  f-om  the 
slot  as  shown  by  figures  20,  21,  and  22  reveals  the  origin  of  the 
propagating  field  as  an  initial  pulse  that  is  formed  from  the 
earliest  and  most  rapidly  increasing  part  of  the  internal  field. 
This  pulse  which  is  obscured  by  the  static  field  in  figure  19  is 
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seen  first  in  figure  20  as  a  spike  that  peels  away  from  the  front 
of  the  static  field  and  arrives  at  the  observation  point  (z  = 
-0.169  m)  with  a  time  delay  of  less  than  1  ns.  However,  the  full 
complexity  of  the  initial  pulse  and  the  reflected  fields  it  gen¬ 
erates  can  be  seen  only  at  observation  points  where  the  static 
field  is  reduced  to  a  level  that  allows  the  use  of  an  expanded 
vertical  scale  as  in  figures  21  to  26.  Figure  21  shows  clearly 
the  basic  sequence  of  events  that  is  repeated  at  every  point  in¬ 
side  the  enclosure.  First  the  initial  pulse  arrives  at  z  -  -0.338 
m  after  a  time  delay  relative  to  the  external  lield  equal  to 
twice  the  time  delay  seen  at  z  =  -0.169.  Approximately  2  to  3  ns 
later,  the  initial  pulse  is  interrupted  by  the  greatly  reduced 
static  field  which  is  interrupted  in  turn  at  t  as  16  ns  by  the 
arrival  of  the  reflected  pulse  generated  by  the  initial  pulse  at 
the  back  wall.  Since  the  reflection  coefficent  at  the  back  wall 
is  +1  ,  the  amplitude  and  polarity  of  this  reflected  pulse  are  at 
first  identical  to  that  seen  in  the  initial  pulse.  However,  at 
later  times  (t  =  18+  ns  and  t  =  22  ns)  pulses  of  opposite  polar¬ 
ity  appear  giving  the  propagating  field  its  characteristic  bipo¬ 
lar  appearance. 

The  presence  of  negative  pulses  is  puzzling  at  first  since 
they  are  obviously  part  of  the  field  reflected  from  the  back  wall 
but  cannot  be  created  by  such  a  reflection.  The  explanation  pro¬ 
vided  by  a  study  of  the  remaining  figures  in  this  series  (espe¬ 
cially  fig.  24,  25,  and  26)  is  that  the  initial  pulse  is  already 
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bipolar  by  the  time  it  reaches  the  back  wall.  This  bipolarity  is 

the  result  of  oblique  reflections  from  the  side  walls  where  the 

boundary  condition  requiring  that  h3x  =  0  at  x  =  +A/2  produces  a 

reflection  coefficient  of  -1.  Thus,  an  initial  pulse  generated 

with  positive  polarity  {same  polarity  as  the  external  field)  at 

_  3 

the  slot  z  =  -6.4  x  10  m  becomes  a  bipolar  pulse  due  to  inter¬ 
actions  with  the  sides  of  the  enclosure  and  is  further  modified 
by  reflections  from  the  front,  back,  and  side  walls  as  it  prop¬ 
agates  in  the  enclosure.  The  negative  sections  or  phases  of  the 
propagating  field  are  labeled  "side-wall"  reflections  in  the  fig¬ 
ures;  however,  most  of  these  are  produced  by  a  series  of  reflec¬ 
tions  involving  the  side  wall  together  with  the  back  and/  or 
front  walls. 

The  earliest  evidence  of  a  negative  phase  in  the  initial 
pulse  is  the  ripple  that  occurs  in  figure  21  just  before  the 
static  field  reaches  its  peak  at  t  a  10  ns.  This  negative  phase 
reappears  at  t  =  18+  ns  after  reflection  from  the  back  wall,  at  t 
=  22  ns  after  reflection  from  the  front  wall,  and  at  later  times 

corresponding  to  the  arrival  of  additional  reflections  from  the 
back  and  front  walls.  The  amplitude  of  the  negative  phase  arriv¬ 
ing  from  the  back  wall  at  t  =  18  ns  is  sufficient  to  cancel  most 


of  the  positive  phase  reflected  from  the  front  wall  which  arrives 
simultaneously  at  this  location. 


The  principal  features  of  the  transient  response  are  seen 
even  more  clearly  in  figure  22  where  z  =  -  0.675  m.  At  this  loca- 
ation,  the  initial  pulse  arrives  after  a  time  delay  equal  to 
twice  the  delay  seen  in  figure  21,  and  its  negative  phase  pro¬ 
duces  a  more  pronounced  ripple  as  it  rides  on  the  rising  portion 
of  the  static  field.  Since  the  location  is  now  closer  to  the  back 
wall  and  farther  from  the  front  wall,  the  back-wall  reflection 
arrives  earlier  and  the  front-wall  reflection  later  than  in 
figure  21.  This  difference  in  arrival  times  effectively  separates 
the  front-wall  reflection  from  the  back-wall  reflection  allowing 
positive  and  negative  phases  of  each  reflection  to  be  seen 
without  interference. 

The  time  delay  of  the  initial  pulse  and  the  separation  of 
the  front-  and  back-wall  reflections  increase  still  more  when  the 
location  moves  to  the  center  of  the  enclosure  (z  =  -1.35  m)  .  As 
shown  in  figure  23,  the  positive  phase  of  the  initial  pulse  ar¬ 
rives  after  a  delay  of  4+  ns,  and  the  first  back-wall  reflection 
arrives  9+  ns  later.  All  subsequent  reflections  of  the  positive 
phase  are  also  spaced  9+  ns  apart.  This  is  exactly  what  one  would 
expect  at  this  location  for  a  pulse  traveling  at  a  constant 
velocity.  Furthermore,  since  the  static  field  is  small,  the 
negative  phase  of  the  initial  pulse  is  clearly  visible  for  the 
first  time  at  t  =  10  ns.  Subsequent  reflections  of  the  negative 
phase  are  indicated  on  the  figure;  however,  these,  unlike  posi- 
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tive  phase  reflections,  do  not  appear  at  completely  regular 
intervals . 

At  the  point  z  =  -2.03  m  midway  between  the  center  of  the 
enclosure  and  the  center  of  the  back  wall  (fig.  24),  the  delay  of 
the  initial  pulse  increases  to  7  +  ns,  and  the  separation  between 
the  initial  pulse  and  its  reflection  from  the  back  wall  is  Bmall 
enough  so  that  its  negative  phase  partially  cancels  the  positive 
phase  of  the  reflection  at  t  =  12  ns. 

Figure  25  shows  the  initial  pulse  arriving  at  the  center  of 
the  back  wall  (z  =  -2.7  m)  with  a  delay  of  9+  ns  and  combining 
with  the  reflected  field  to  give  a  net  field  with  the  same  form 
as  the  incident  field,  The  initial  pulse  is  interrupted  at  t  =  27 
ns  by  the  return  of  its  positive  phase  after  reflection  from  the 
front  wall.  Succeeding  reflections  from  the  front  wall  arrive  at 
intervals  of  18+  ns  corresponding  to  the  time  required  for  the 
pulse  to  travel  from  the  back  wall  to  the  front  wall  and  return. 
The  18  ns  interval  between  the  arrival  of  the  initial  pulse  and 
the  return  of  its  positive  phase,  combined  with  the  absence  of 
the  static  field,  provides  the  clearest  picture  of  the  initial 
pulse  to  be  seen  in  these  figures.  It  shows  that  the  pulse  in¬ 
cludes  not  just  a  single  cycle  of  positive  and  negative  phases 
but  also  a  second  cycle  with  reduced  magnitudes  beginning  at  t  = 
20  ns.  Evidently,  the  initial  pulse  consists  of  a  train  .  f  such 
cycles  generated  by  side-wall  reflections.  The  folding  of  this 
train  back  onto  itself  at  each  reflection  from  the  front  and  back 
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walls  accounts  for  the  increasing  complexity  with  time  of  the 
transient  waveform. 

Since  the  reflection  coefficient  of  the  back  wall  is  +1,  we 
should  expect  the  initial  pulse  to  be  doubled  in  magnitude  when 
it  combines  with  its  reflection  at  z  =  -2.7  m  in  figure  25.  One 
way  to  verify  this  result  is  to  recompute  the  field  at  z  =  -2.7  m 
for  the  same  enclosure  with  the  back  wall  removed.  This  has  been 
done  in  figure  26  by  setting  C  =  5.4  m,  that  is,  by  analytically 
moving  the  the  back  wall  to  z  =  -5.4  m.  Comparing  figure  25  to 
26,  we  see  that  the  peak  magnitude  of  the  initial  field  in  figure 
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6. 


COMPARISON  BETWEEN  COMPUTED  AND  MEASURED  FIELDS  INSIDE 


SHIELDED  ENCLOSURES 

Although  rectangular  enclosures  have  been  used  for  at  least 
sixty  years  as  electromagnetic  shields,  there  is  surprisingly 
little  experimental  data  in  the  avaliable  literature  that  is 
directly  relevant  to  this  type  of  shield.  Most  published  data 
relates  to  shields  with  simpler  geometry,  such  as  the  plane  sheet 
and  the  sphere.  This  is  probably  due  to  a  desire  to  avoid  the 
complexity  of  the  rectangular  geometry  by  using  more  amenable 
structures  that  might  be  expected  to  approximate  the  rectangular 
enclosure  for  certain  ranges  of  frequency,  for  example,  the  sheet 
for  very  high  frequencies  and  the  sphere  for  very  low  frequen¬ 
cies.  Moreover,  those  reports  that  do  provide  data  on  rectangular 
enclosures  frequently  fail  to  disclose  all  the  relevant  informa¬ 
tion.  Examples  of  commonly  omitted  information  are  the  enclosure 
dimensions,  the  location  of  the  internal  field  sensor,  and  the 
locations  and  orientation  of  slots  and  seams. 

In  this  section,  we  reproduce  two  of  the  most  complete  sets 
of  experimental  shielding  data  for  rectangular  enclosures  and 
compare  this  data  to  calculations  based  on  the  theoretical  expes- 
sions  we  have  developed  in  the  preceding  sections.  Both  sets  are 
frequency  domain  data  and,  as  is  the  usual  practice  for  measure¬ 
ments  of  this  type,  are  expressed  in  terms  of  the  shielding 


effectiveness  (or  shielding  efficiency)  SE  of  the  enclosure.  SE 
is  defined  as  follows: 


SE 


20 


(dB) 


(6.1) 


where  is  the  pth  component  of  the  magnetic  field  measured  at 

a  specified  location  inside  the  enclosure  generated  by  a  specifi¬ 
ed  source  outside  the  enclosure  and  H°  is  the  same  component  gen¬ 
erated  by  the  same  source  at  the  same  location  when  the  enclosure 
walls  and  ceiling  have  been  removed.  That  is,  H  is  the  shielded 

p 

field  at  a  location,  H°  is  the  unshielded  field  at  the  same  loca- 

P 

tion  due  to  the  same  source,  and  the  logarithmic  difference  be¬ 
tween  the  magnitudes  of  these  quantities  in  dB  is  defined  as  the 
shielding  effectiveness  of  the  enclosure  at  that  location.  This 
definition  is  intended  to  characterize  the  shielding  properties 
of  the  enclosure  in  a  way  that  is  independent  of  the  source  ex¬ 
cept  for  its  frequency.  If  SE  can  be  measured  in  a  way  that  is 
independent  of  the  source,  then  it  can  be  used  to  predict  the 
shielding  performance  of  the  enclosure  against  the  magnetic  field 
from  any  source,  provided  that  the  frequency  of  source  is  known. 
However,  in  general,  it  is  not  possible  to  measure  SE  at  any 
point  inside  an  enclosure  in  a  way  that  is  independent  of  the 
source.  This  is  due  to  the  fact  that  two  different  sources  S^  and 
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Sg  will  produce  unshielded  fields,  say  and  H”“ ,  that  vary  in 
different  ways  throughout  the  volume  such  that  in  general  at  any 
location 


IHpl/IMp'l  '  |Hp|/|H“2|  (6.2) 

1  2 

where  H  and  H  are  the  shielded  fields  corresponding  to  and 

P  P  ^ 

Sg  respectively.  Nevertheless,  there  are  classes  of  sources  for 


which 
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(6.3) 


for  any  sources  Sj  and  belonging  to  such  a  class.  In  this 
case,  when  SE  is  obtained  from  (6.1)  using  measurements  made  with 
one  source  belonging  to  this  class,  it  characterizes  the  shield¬ 
ing  of  the  enclosure  with  respect  to  all  sources  in  that  class. 
That  is,  once  SE  is  determined  in  this  way,  it  can  be  used  to 
predict  the  shielding  properties  of  the  enclosure  against  magnet¬ 
ic  fields  from  any  source  belonging  to  this  class.  An  example  of 
a  class  of  sources  for  which  (6.3)  holds  is  the  class  of  all 


sources 


that  generate  spatially  uniform  unshielded  fields 


throughout  the  region  of  interest.  For  any  two  sources  and  S£ 


belonging  to  this  class 

Hp2 


(6.4) 


at  each  point  where  AQ(jw)  may  be  a  constant  or  a  function  of 
frequency  but  is  independent  of  the  spatial  variables.  In  this 
case,  the  shielded  fields  will  obey  the  same  relationship: 

H1  =  A  (jw)  H2  (6.5) 

p  o w  p 

at  each  point  and  (6.3)  will  be  satisfied.  Consequently,  SE  as 


defined  by  (6.1)  completely  characterizes  the  shielding  proper¬ 
ties  of  the  enclosure  with  respect  to  the  class  of  sources  that 
expose  the  enclosure  to  spatially  uniform  fields. 

We  can  use  our  expression  for  the  field  inside  a  continuous 
enclosure  exposed  to  uniform  fields  on  one  wall  to  compute  the 
theoretical  shielding  effectiveness  of  these  structures  against 
the  class  of  uniform  sources  by  substituting  H3^ ( x , y , z ; ju )  from 

(4.6)  for  the  shielded  field  H  and  HlS(jw)  for  the  unshielded 

P  P 

field  H°  in  (6.1): 

P 


| H 3  ( x , y , z ; jw  )  | 

SE(x,y,z;w)  =  -20  log  - £ - 
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(6.6) 


(When  more  than  one  wall  is  exposed  to  the  field  the  contributios 
from  all  must  be  added  at  the  specified  location  to  determine  the 
total  unshielded  field  H3^( x , y , z ; jw ) ) .  Calculations  using  (6.6) 
can  then  be  compared  with  measurements  based  on  (6.1)  to  deter¬ 
mine  whether  or  not  the  theory  is  supported  by  experiment.  It 
should  be  noted  from  (6.6)  that  SE  is  a  function  of  the  spatial 
variables  even  though  the  unshielded  field  is  not. 


Continous  Copper  Enclosure  with  Soldered  Joints 


In  1964,  J.  B.  Hays  published  the  results  of  shield- 
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ing  measurements  that  had  been  carried  out  by  E.  D.  Sunde  on  an 
8  x  8  x  8  ft.  cubical  structure  constructed  from  10  mil  copper 
sheets  soldered  together  along  their  edges.  The  experimental  set¬ 
up  used  in  these  measurements  is  represented  schematically  in 
part  A)  of  figure  27.  It  consists  of  a  signal  generator  driving  a 
circular  coil  80  ft.  in  diameter  surrounding  the  copper  cube 
which  contains  a  smaller  coil  used  as  a  magnetic  field  detector. 
The  purpose  of  the  large  coil  and  generator  is  to  produce  an 
unshielded  magnetic  field  with  a  single  spatially  uniform  compon- 
nent  H°  over  the  volume  occupied  by  the  enclosure.  This  allows 
easy  determination  of  H°  .which  can  then  be  combined  with  the 

y 

shielded  field  measured  by  the  loop  detector  to  give  SE  using 
(6.1).  With  this  setup,  Sunde  was  able  to  measure  the  shielding 
effectiveness  of  the  copper  cube  over  a  frequency  range  from  less 
than  100  Hz  to  40  kHz.  His  data  (reproduced  from  curve  d  of  fig. 
1  in  reference  18)  is  shown  as  the  broken  line  in  part  B)  of  fig¬ 
ure  2  7  . 

The  solid  curve  in  part  B)  is  the  computed  shielding 
effectiveness  obtained  with  the  aid  of  H3^(y,z;jw)  from  equation 
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it  was  assumed  that  the  detector  was  located  at  the  center  of  the 
cube  which  in  our  standard  coordinate  system  is  given  by  x  =  y  = 
0,  z  =  -1.219  m.  (The  actual  location  was  not  mentioned  by  Hays, 
but  it  appears  highly  likely  that  a  desire  to  preserve  symmetry 
would  have  dictated  this  choice).  Since  the  total  field  at  the 
center  of  a  uniformly  illuminated  cube  is  the  sum  of  equal  con¬ 
tributions  from  fields  entering  through  each  of  the  cube’s  four 
vertical  walls,  the  theoretical  expression  for  the  shielded  field 
at  this  location  is  4  H3  ( 0 , - 1 . 2 1 9 ; ju )  and  the  computed  shielding 
effectiveness  is  obtained  by  evaluating 


SE ( 0 , 0 , - 1 . 2 1 9 ; w )  =  -20  log 


as  a  function  of  f  =  u/2n  for 


4  |H3  (0,-1.219;  ju)  | 

|HlJ<j«)  | 


(6.7) 


A  =  B  =  C  =  2.438  m,  d  =  2.54x10  m 
o2  =  5.8x10^  mhos/m,  p 2  =  12.6x10  ^h/m. 

Comparison  of  the  two  curves  shows  good  agreement  at  the 
higher  frequencies.  There  is,  however,  a  divergence  below  4  kHz 
that  causes  the  measured  values  of  SE  to  fall  below  the  computed 
values.  This  divergence  may  be  due  to  the  fact  that  condition 
(2.8)  is  not  satisfied  at  these  frequencies  for  10  mil  copper 
sheets.  That  is,  the  thickness  of  the  sheet  is  not  larger  than 
the  skin  depth  of  the  field  in  the  sheet  at  these  frequencies  as 
assumed  by  the  theory.  This  means  that  multiple  reflections  can 
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occur  inside  the  sheet  which  may  result  in  a  larger  internal 
field  and  a  smaller  shielding  effectiveness  than  can  be  accounted 
for  by  the  theory  in  its  present  form. 

6.2  Slotted  Steel  Enclosure 


The  effect  of  slot  length  a  on  the  shielding  effec¬ 
tiveness  of  a  cubical  steel  enclosure  was  investigated  by  E.  M. 

20 

Honig,  Jr.  using  the  experimental  setup  shown  in  part  (a)  of 

figure  28.  Unlike  the  arrangement  employed  by  Sunde  for  the 
continuous  enclosure,  this  setup  is  not  designed  to  generate  a 
spatially  uniform  source  field.  Instead,  a  Bmall  loop  source  is 
used  to  generate  a  strongly  nonuniform  field  throughout  the  vol¬ 
ume  occupied  by  the  enclosure.  The  advantages  of  Honig* s  setup 
are  that  it  can  be  used  at  higher  frequencies  and  generally  pro¬ 
duces  a  greater  dynamic  range  for  measurements  on  enclosures  of 
this  size  than  does  Sunde’s.  Its  disadvantage  is  that  measure¬ 
ments  of  SE  made  in  this  way  using  (6.1)  depend  on  the  location 
of  the  source.  That  is,  the  same  source  located  at  a  different 
position  with  respect  to  the  enclosure  and  the  detector  would 
produce  a  different  value  of  SE.  Similarly,  a  larger  or  smaller 
source  at  the  same  position  might  produce  a  different  value  of 
SE.  To  mitigate  this  limitation,  Honig  has  chosen  the  setup 
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called  for  by  two  widely  used  standards  for  shielding  measure- 
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ments:  MIL-STD  285  and  IEEE  299  which  specify  that  the 

source  and  the  detector  shall  be  identical  circular  loop  antennas 
12  in.  in  diameter  lying  in  the  same  plane  (coplanar 
arrangement)  with  centers  18  in.  from  the  enclosure  wall.  With 
this  set-up,  Honig’s  measurements  can  be  compared  directly  to 
other  measurements  made  using  these  same  standards,  even  though 
they  would  not  in  general  be  comparable  to  measurements  made  on 
the  same  enclosure  using  a  different  source  such  as  Sunde’s  large 
loop.  In  accordance  with  these  standards,  Honig  orients  the  plane 
of  his  loops  perpendicular  to  the  long  axis  of  the  slot.  This 
orientation  produces  the  most  efficient  coupling  through  slots  of 
any  size  and  insures  that  his  measured  values  of  shielding 
effectiveness  will  be  the  smallest  possible  for  this  setup.  It 
allows  him  to  measure  the  shielding  effectiveness  of  slots 
ranging  from  1  to  30.5  cm  in  length  at  frequencies  of  10  kHz,  40 
kHz,  200  kHz,  and  1  MHz  for  each  slot.  The  procedure  is  to  remove 
the  enclosure  wall  containing  the  slot  and  measure  the  unshielded 
field  H°  at  the  detector  loop  due  to  the  source  loop  and  then  to 

replace  the  wall  and  measure  the  shielded  field  H  at  the 

r  x 

2 1 

Dept  of  Defense,  MIL-STD-285,  Method  of  Attenuation 
Measurements  for  Enclosures,  Electromagnetic  Shielding,  for 
Electronic  Test  Purposes  (25  June  1956). 
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detector  without  changing  the  arrangement  of  the  loops  from  that 
shown  in  the  figure.  The  latter  step  is  then  repeated  as  the  slot 
length  is  increased  in  specified  increments.  The  measured 
shielding  effectivness  for  each  length  is  determined  by  (6.1). 

Honig’s  data  from  figure  8  in  his  report  is  plotted  versus 

slot  length  in  part  (b)  of  figure  28  along  with  a  data  point  ob- 

23 

tained  by  Axford  et  al.  who  measured  the  shielding  effecetive- 
ness  of  a  cubical  steel  enclosure  with  a  50  cm  slot  using  a  very 
similar  setup.  The  data  shows  a  strong  general  trend  in  the  di¬ 
rection  of  decreasing  shielding  effectiveness  with  increasing 
slot  length.  However,  there  is  a  notable  subset  of  the  data  (a  = 
7.5  ,  10,  and  13.5  cm)  where  this  trend  is  reversed.  The  plot 
also  shows  that,  with  the  exception  of  one  point  (a  =  7.5  cm), 
the  data  for  each  slot  is  closely  grouped  indicating  that  the 
shielding  effectiveness  is  virtually  independent  of  the  source 
frequencies  used  in  these  measurements. 

Although  Honig’s  loop  source  generates  a  nonuniform  un¬ 
shielded  field  H°(x,y,z;jw)  throughout  most  of  the  volume  occu¬ 
pied  by  the  slotted  enclosure,  it  will  produce  a  nearly  uniform 


R.  Axford,  R.  McCormack,  and  R.  Mittra,  Evaluation  of  the 
Applicability  of  Standard  CW  EMI/RFI  Shielding  Effectiveness  Test 
Techniques  to  Assessment  of  EMP  Hardness  of  Tactical  Shelters, 
Construction  Engineering  Research  Laboratories,  CERL-TM-M-307 
Urbana ,  IL  (March  1982). 
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field  over  a  slot  positioned  like  the  one  shown  in  the  figure 
provided  the  slot  is  not  too  long.  For  such  a  slot,  we  can  calc- 
culate  the  field  H3x { x , y , z ; jw )  at  any  point  inside  the  enclosure 
in  terms  of  the  tangential  compnent  of  the  unshielded  loop  field 
H°(0,0,0;jw)  incident  on  the  slot  by  setting  H°(0,0,0;jw)  = 
Hl^(jw)  in  equation  (4.24).  With  this  substitution,  the  shielded 
field  at  the  center  of  the  detector  loop  (0,0, -0.49  m)  due  to  the 
source  loop  is  H3x ( 0 , 0 , -0 . 49  m;ju),  and  the  theoretical  expres¬ 
sion  for  the  shielding  effectivness  that  would  be  measured  at 

that  location  using  Honig’s  set-up  is 

'  |H3  (0,0,-.49m;jw)  |  ' 

SE ( 0 , 0 , - . 49m ; u )  =  -20  log  - - -  (6.8) 

|Hjc  (  ©  *  ©  »  ~  •  49m ;  ju> )  | 

where  H°(0,0,-0.49  m;ju)  is  the  unshielded  source  field  at  the 
center  of  the  detector  loop.  Since  H3x ( 0 , 0  , -0 . 49  m;jw)  is  comput¬ 
ed  in  terms  of  H° ( 0 , 0 , 0 , ; jw ) ,  it  is  also  convenient  to  write  the 
field  H°(0,0,-0.49  m;jw)  in  terms  of  H° ( 0 , 0 , 0 ; jw  )  .  This  is  easi¬ 
ly  done  using  the  fact  that  points  (0,0,0)  and  (0,0, -0.49  m)  are 
both  in  the  near  field  of  the  source  loop  where  the  magnitude  of 
the  field  is  inversely  proportional  to  the  cube  of  the  distance 
to  the  center  of  the  loop* .  In  this  case,  since  (0,0, -0.49  m)  is 
twice  as  far  from  the  center  of  the  source  loop  as  (0,0,0),  the 

*R.  L.  Monroe,  A  Theory  of  Electromagnetic  Shielding  with 
Applications  to  MIL-STD  285,  IEEE-299,  and  EMP  Simulation, 

Harry  Diamond  Laboratories,  HDL-CR-8 5-0 5 2- 1 ,  Adelphi ,  MD 
( February  1985 ) . 
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magnitude  of  Hx(0,0,-0.49  m;jw)  will  be  equal  to  1/8  that  of 
H° ( 0 , 0 , 0 ; jw ) .  Thus  (6.8)  can  be  rewritten 


SE ( 0 , 0 , - . 49m ; w )  =  -  20  log 


8  |H3  (0,0t-.49m;  jtt)  | 
|H°(0,0,0;jw)  | 


(6.9) 


Since  H3x ( 0 , 0 , -0 . 49  m;ju)  is  directly  proportional  to 
H° ( 0 , 0 , 0  ;  jw )  ,  this  term  drops  out  when  (6.9)  is  evaluated;  and  SE 
is  seen  to  depend  on  the  source  only  through  the  factor  8. 

The  calculated  shielding  effectivness  for  Honig’s  setup  is 
denoted  by  the  solid  triangles  at  a  =  1.56,  3.13,  6.25,  12.5,  25, 
and  50  cm  in  part  (b)  of  figure  28.  This  plot  was  obtained  by 
evaluating  (6.9)  for  each  a  (expressed  in  meters)  at  a  frequency 
of  10  kHz.  Only  a  single  frequency  was  required  because  previous 
calculations  (fig.  11  and  12)  based  on  equation  (4.24)  demonstra¬ 
ted  that  H3x (x , y , z ; jw )  is  independent  of  frequency  for  frequen¬ 
cies  below  the  cutoff  frequency  of  the  enclosure.  The  enclosure 
parameters  used  in  these  calculations  are 
A  =  B  =  C  =  2.7  m  ,  d  =  6.4xl0-3  m 
o2  =  4x10^  mhos/m,  p2  =  12.6x10'^  h 
and  the  slot  inductance  for  each  slot  length  a  (meters)  was 

computed  with  equation  (4,29): 

-  8 

L  =  6.4  x  10  a  (henries) 


14 

which  was  derived  from  figure  7  of  Monroe  .  Since  Honig  did  not 

specify  the  location  of  the  slots,  it  was  assuned  that  they  are 

centered  at  the  center  of  the  wall  where  x  =  y  =0. 

c  c 

Comparison  of  the  measured  and  computed  shielding  effec¬ 
tiveness  curves  shows  good  agreement  for  slots  less  than  6  cm  in 
length.  A  divergence  occurs  at  a  =  7.5  cm,  where  the  measured 
shielding  effectiveness  is  larger  than  the  preceding  measurments 
at  5.5  and  6.5  cm  rather  than  smaller  as  shown  by  the  computed 
curve.  Similar  divergences  occur  at  a  =  10  and  13.5  cm.  These 
data  points  run  counter  to  the  general  trend  of  both  measured  and 
computed  curves  and  are  difficult  to  explain  on  physical  grounds 
since  they  would  require  a  type  of  antiresonant  behavior  that 
appears  to  be  ruled  out  at  the  frquencies  used  in  the  measure¬ 
ments.  Consequently,  the  validity  of  these  three  points  appears 
to  be  open  to  question.  If  they  were  eliminated,  the  measured 
curve  like  the  computed  curve  would  show  SE  as  a  monotonically 
decreasing  function  of  a,  and  the  overall  agreement  would  be 
noticeably  improved.  For  very  small  slots  (<  1  cm),  SE  exceeds 
the  dynamic  range  of  Honig’ s  measurements  (100  to  125  dB ) ,  and, 
therefore,  valid  comparisons  cannot  be  made.  However,  the  figure 
suggests  that  SE  will  continue  to  increase  monotonically  as  slot 

*^R.  L.  Monroe,  EMP  Shielding  Effectiveness  and  MIL-STD  285, 
Harry  Diamond  Laboratories,  HDL-TR-1336,  Adelphi,  MD.  (July 
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length  decreases,  until  it  approches  the  shielding  effectiveness 


of  an  enclosure  with  no  slot,  that  is,  until  it  approaches  the 


shielding  effectiveness  of  a  continuous  steel  enclosure.  To  calc¬ 


ulate  this  limit,  it  would  be  necessary  to  use  a  modified  version 


of  (6.9)  where  h3^ ( 0 , 0 , -0 . 49  m;ju)  is  computed  with  an  expression 


similar  to  (4.6). 


7. 


CONCLUSIONS 


In  the  preceding  pages,  we  have  used  a  theory  of  electro¬ 
magnetic  shielding  based  on  impedance  boundary  conditions  to  ob¬ 
tain  general  expressions  for  the  fields  at  any  point  inside  slot¬ 
ted  and  continuous  rectangular  enclosures  exposed  on  one  side  to 
arbitrary  external  source  fields.  These  expressions  consist  of 
infinite  series  summed  over  the  TE  waveguide  modes  where  each 
term  is  the  product  of  one  or  more  mode  functions  and  a  Fourier 
coefficient.  In  the  case  of  the  slotted  enclosure,  the  Fourier 
coefficients  are  determined  by  the  spatial  distribution  of  the 
magnetic  field  incident  on  the  inside  surface  of  the  wall  due  to 
the  slot,  and  for  the  continuous  enclosure  they  are  determined  by 
the  spatial  distribution  of  the  magnetic  field  incident  on  the 
outside  surface  of  the  wall.  These  coefficients  were  evaluated  in 
closed  form  for  the  case  where  the  outside  surface  of  the  slot  is 
exposed  to  a  spatially  uniform  magnetic  field  and  for  the  case 
where  the  entire  outside  surface  of  the  continuous  enclosure  is 
exposed  to  a  spatially  uniform  magnetic  field.  Frequency  domain 
calculations  based  on  these  expressions  showed  that  the  series 
converge  rapidly.  This  allowed  easy  plotting  of  the  internal 
magnetic  fields  as  functions  of  frequency  with  the  aid  of  a 


microcomputer.  The  plots  show  that  at  frequencies  below  the  cut¬ 
off  frequency  the  field  inside  a  typical  slotted  enclosure  is 
independent  of  freauency  and  the  field  inside  a  continuous 
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enclosure  is  a  monotonies lly  decreasing  function  of  frequency.  At 
frequencies  above  cutoff,  the  fields  in  both  slotted  and  contin¬ 
uous  enclosures  show  very  complex  behavior  associated  with  cav¬ 
ity  resonances.  The  spatial  variation  of  the  internal  fields  is 
determined  by  the  boundary  conditions  at  the  inside  surfaces  of 
the  walls.  The  field  reaches  its  maximum  at  the  exposed  wall  and 
decreases  monotonically  as  one  moves  to  the  opposite  wall. 

The  frequency  domain  expressions  were  replaced  by  series 
approximations  and  transformed  into  the  time  domain  using  a  term 
by  term  application  of  the  inverse  LaPlace  transform.  The  trans¬ 
formed  series  were  found  to  converge,  and  they  were  used  to  com¬ 
pute  and  plot  the  transient  response  of  the  internal  magnetic 
fields  when  one  wall  of  the  enclosure  is  exposed  to  time-varying 
external  fields  in  the  form  of  decaying  exponentials,  damped 
sinusoids,  or  rationalized  exponentials.  These  plots  showed  quite 
different  responses  for  the  continuous  and  slotted  enclosure  in 
agreement  with  the  results  of  Monroe  *  who  computed  the  transient 
fields  at  one  point  on  the  inside  surface  of  the  exposed  wall 
using  one-term  approximations.  It  was  found  that  the  character¬ 
istics  of  the  field  inside  a  continuous  enclosure  are  determined 
to  a  great  extent  by  the  diffusion  time  factor  r  (eq  5.4)  and  by 

^R.  L.  Monroe,  A  Theory  of  Electromagnetic  Shielding  with 
Applications  to  MIL-STD  285,  IEEE-299,  and  EMP  Simulation, 

Harry  Diamond  Laboratories,  HDL-CR-85 -052- 1 ,  Adelphi,  MD 


(  February  1985)  . 
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the  duration  of  the  incident  field.  The  former  determines  the 
shortest,  that  is,  fastest,  rise  time  possible  by  the  internal 
field,  and  the  latter  determines  its  largest  value.  Plots  of  the 
fields  inside  slotted  enclosures  showed  a  transient  response  con¬ 
sisting  of  the  sum  of  stationary  and  propagating  fields.  The  sta¬ 
tionary  field  dominates  close  to  the  slot  where  it  forms  a  reduc¬ 
ed  replica  of  the  external  field.  However,  its  peak  field  de¬ 
creases  rapidly  with  distance  from  the  slot  so  that  at  most  loca¬ 
tions  in  the  interior  it  is  negligible  compared  to  the  propagat¬ 
ing  field.  The  latter  travels  at  approximately  the  speed  of  light 
and  undergoes  repeated  reflections  from  the  front,  back,  and 
sides  of  the  enclosure.  These  reflections  produce  a  waveform  that 
becomes  increasingly  complex  with  time.  The  separation  of  the 
internal  field  into  these  two  constituents  corresponds  to  a  sep¬ 
aration  of  the  incident  field  into  the  sum  of  frequencies  lying 
above  and  below  the  cutoff  frequency  with  the  propagating  field 
formed  from  frequencies  above  cutoff  and  the  stationary  field 
formed  from  those  below.  The  peak  value  of  the  stationary  field 
depends  principally  on  the  length  of  the  slot,  while  the  peak  of 
the  propagating  field  depends  on  the  rise  time  of  the  external 
field.  That  is,  longer  slots  produce  larger  stationary  fields  and 
faster  rise  times  produce  larger  propagating  fields. 


* 
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Finally,  oalcu’ations  based  on  the  frequency  domain  expres¬ 
sions  were  used  to  compare  the  theoretical  results  with  two  sets 

18  19 

of  experimental  data:  one  set  by  E.  D.  Sunde  ’  for  a  contin¬ 
uous  cubical  enclosure  formed  from  soldered  copper  sheets  and  the 

20 

other  by  E.  M.  Honig  for  a  cubical  steel  enclosure  with  one 
wall  penetrated  by  single  slot  with  a  variable  length.  Since  the 
data  was  presented  in  terms  of  the  measured  shielding  effective¬ 
ness  SE  of  each  enclosure,  it  was  necessary  to  construct  theoret¬ 
ical  expressions  for  SE  corresponding  to  the  experimental  setups 
of  both  Sunde  and  Honig  using  H3  (y,z;ju)  from  equation  (4.6)  and 

y 

H3^ ( x , y , z ; ju )  from  equation  (4.24).  Computed  values  of  SE  were 
then  plotted  together  with  the  experimental  values  for  easy 
comparison.  In  the  case  of  the  copper  enclosure,  good  agreement 
is  found  for  frequencies  above  4  kHz,  but  a  significant  diverg¬ 
ence  occurs  at  extremely  low  frequencies.  For  the  slotted  enclo¬ 
sure,  there  is  good  agreement  for  slots  less  than  6  cm  in  length, 
but  a  divergence  is  seen  for  longer  slots  where  some  of  the 
experimental  data  points  show  longer  slots  providing  greater 
shielding  than  some  shorter  slots.  Since  good  agreement  is  seen 

B.  Hays,  IEEE  Spectrum,  Vol  1  ,  May  (1964). 
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E.  D.  Sunde,  Switching  Center  Shielding  Against  Atmospheric 
Induction,  Bell  Telephone  Laboratories  (unpublished  memorandum). 
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where  the  theory  is  expected  to  be  most  accurate  (high  frequen¬ 
cies  for  continuous  enclosures  and  short  slot  lengths  for  slotted 
enclosures),  the  overall  comparison  provides  a  satisfactory  veri¬ 
fication  of  the  theory  within  its  present  limits. 

These  results  establish  a  theoretical  basis  for  the  princi¬ 
pal  features  of  electromagnetic  shields  formed  from  rectangular 
enclosures.  However,  they  do  not  exhaust  the  subject.  There  re¬ 
main  many  areas  open  for  further  investigation  relating  to  both 
the  general  shielding  theory  and  its  application  to  practical 
shielded  enclosures. 

One  area  for  possible  future  development  is  an  extension 
of  the  general  shielding  theory  to  include  thin-walled  structures 
where  conditions  (2.8)  and  (2.20)  ar  not  satisfied.  Such  an  ex¬ 
tension  might  eliminate  the  divergence  between  measured  and  com¬ 
puted  SE  for  the  thin-walled  copper  enclosure  at  extremely  low 
frequencies  and  clear  the  way  for  accurate  calculations  at  still 
lower  frequencies.  It  might  also  close  the  gap  between  theory 
and  measurement  for  enclosures  with  long  slots  and  lead  to  cor¬ 
responding  results  for  an  enclosure  with  a  resonant  slot. 
(Although  the  effect  of  a  resonant  slot  can  be  partially  account¬ 
ed  for  simply  by  replacing  q 2  or  q 2  in  (3.35)  with  a  general 


r.  %  _  *  .  *-r.  «r.  »_  itl  *_  r.  «r\.  . 


expression  for  the  slot  impedance  ,  this  would  not  change  the 
fact  that  condition  (2.20)  would  not  be  satisfied,  and  conse¬ 
quently  the  accuracy  of  the  resulting  expression  could  not  be 
guaranteed).  To  make  this  extension  with  assurance,  it  would  be 
necessary  to  replace  the  simple  impedance  boundary  conditions 
(2.2)  and  (2.10)  with  more  general  forms  that  account  for  multi¬ 
ple  reflections  inside  the  walls  and  slots  of  the  enclosure.  If 
these  boundary  conditions  can  be  developed,  then  there  appears  to 
be  no  reason  why  the  remainder  of  the  extension  cannot  be  carried 
out  in  nearly  the  same  manner  as  the  original.  The  extension 
would,  of  course,  be  more  complicated  than  the  original,  but  its 
likely  advantage  would  compensate  for  the  added  complexity 

An  extension  of  the  theory  to  enclosures  with  slots  filled 
with  lossy  dielectrics  could  also  be  very  useful.  This  extension 
would  include  both  the  inductive  (air  filled)  and  resistive  slots 
as  special  cases  and  would  be  a  natural  model  for  the  radio  fre¬ 
quency  gasket-  a  device  that  is  often  used  in  practical  enclo¬ 
sures  to  reduce  the  magnitude  of  unwanted  fields  in  open  seams. 
It  could  also  be  used  to  model  the  effect  of  paints  and  metal 
oxides  that  tend  to  build  up  in  any  unintended  opening.  These 
results  could  be  incorporated  directly  into  the  theory  in  its 
present  form  by  regarding  the  slot  as  a  rectangular  waveguide 
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filled  with  a  lossy  dielectric  and  adapting  the  existing  theory 
of  these  structures  to  the  shielding  problem. 

A  reader  who  has  followed  the  development  to  this  point 
will  no  doubt  be  able  to  suggest  many  other  areas  for  further 
research  on  this  interesting  subject. 
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APPENDIX  A-PARTIAL  FRACTION  EXPANSIONS  OF  Rj(z;ju)  AND  R£m(z;jw) 


Our  objective  here  is  to  obtain  partial  fraction  expansions 

of  Rgfzjjw)  and  R^^zjjw)  which  are  defined  by  equation  (5.36)  on 

page  106.  This  procedure  can  be  shortened  by  noting  that  both  of 

these  quantities  can  be  written  ir  the  form: 

2  ip  {  cosh  [  $  (  C+d  +  z  ) /C  ] 


R(z;$)  = 


( A- 1  ) 


if  <  coshU  )  +  sinh<<) 

where  $  is  a  complex  variable  and  if  is  a  real  parameter.  To 
recover  R^zjjw)  and  R^tzjjw)  from  (A-l),  it  is  only  necessary 


to  replace  (  by 


S  =  (C/CQ) (w‘ )2  +  ( ju)2 


*  by 


L  C 

y  o 

rj  C 
'  o 


(A-2) 


( A-  3  ) 


and  wl  by  wn  and  wnm  respectively, 
c  c  c  r 

Since  R(z;$)  is  a  single  valued  function  of  $  whose  only 
singularities  are  simple  poles,  it  satisfies  the  definition  of  a 
meromorphic  function  and  can  be  expanded  as  a  series  of  partial 


fractions  in  the  form 


R  (  z  5  $  )  =  R  (  z  ;  0  )  +  Wk  - - -  +  — - 

k  ?  "  ?k  ?k  • 


( A-4  ) 


AE.  C.  Titchmarsh,  The  Theory  of  Functions,  2d  Ed.,  Oxford 
University  Press  (1939). 
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where  <k  is  the  kth  pole  of  R(z;$)»  is  its  residue,  the  sub  is 


taken  over  all  poles,  and 
R(z;0)  =  2  f-/(f+l) 

Now  the  poles  of  R(z;$)  are  the  nonzero  roots  of 
F($)  =  t  {  cosh($)  +  sinh($) 


(A-5) 


( A-6 ) 


and  it  can  easily  be  shown  that  these  form  a  doubly  infinite 


sequence  of  paired  imaginary  numbers 


«k  *  ♦j-’k 

where  p,  satisfy 


k  =  1,2,3, 


f  p  cos(p)  +  sin(p)  =  0  . 


(A-7) 


( A-8 ) 


Thus  (A-4)  becomes 


R( z  ;  $  )  =  R(z;0) 


oo 

♦  Kbr 

k  =  l  * 


J Pu  JPl 


( A-9  ) 


5 «;  [— * - - 

-  -  <  +  J Pl-  JPl 


where  are  the  residues  for  the  +  jpk  poles  and  are  the 

residues  for  the  -jpk  poles.  Evaluating  the  residues,  we  find 


Wk ( z )  =  -  Wk ( z )  = 


2*jpk  cos[pk(C+d+z)/C] 
?[cos(p,  )-p,  sin(p,  )  ]  +  cos (p.  ) 


( A- 1 0 ) 


Using  the  first  equation  in  (A-10),  we  can  combine  the  two  sums 
in  (A. 9)  as  follows 


go 

:(z;$)  =  R(z;0)  +  ^w*(z)  - - - 

k=l  «  “ 


$  +  JPl 


( A- 1 1 ) 


W*(z) 


R  (  z  ,  0  )  +  2  $2  \  - 2—^ - g- 

k=l  ^k[<  -  (jpk)  ] 


The  latter  can  be  rewritten  using  (A-2),  (A-3),  and  (A-10),  as 


f ollows 


R( 2 ; jw )  = 


where 


,  2  21  V  Wk(z> 

+  2  (w  )  +  (  jw  )  /  - 7 — 2 - 2 - 

°  ktl  <«c>  +  <V  + 


Wk(z)  =  W*(z)/(jpk) 
wk  =  pkCo/C 

and  the  partial  fraction  expansions  of  and  R£  (z»jw^  can 

„  ,  .  t  ...  n  j  nm 

be  obtained  from  the  preceding  by  replacing  wc  with  wc  and  uc  . 
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